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Abstract. Astrophysical fluids are turbulent a fact which 
changes the dynamics of many key processes, including mag- 
netic reconnection. Fast reconnection of magnetic field in 
turbulent fluids allows the field to change its topology and 
connections. As a result, the traditional concept of magnetic 
fields being frozen into the plasma is no longer applicable. 
Plasma associated with a given magnetic field line at one in- 
stant is distributed along a different set of magnetic field lines 
at the next instant. This diffusion of plasmas and magnetic 
field is enabled by reconnection and therefore is termed "re- 
connection diffusion". The astrophysical implications of this 
concept include heat transfer in plasmas, advection of heavy 
elements in interstellar medium, magnetic field generation 
etc. However, the most dramatic implications of the concept 
are related to the star formation process. The reason is that 
magnetic fields are dynamically important for most of the 
stages of star formation. The existing theory of star forma- 
tion has been developed ignoring the possibility of reconnec- 
tion diffusion. Instead, it appeals to the decoupling of mass 
and magnetic field arising from neutrals drifting in respect to 
ions entrained on magnetic field lines, i.e. through the pro- 
cess that is termed "ambipolar diffusion". The predictions of 
ambipolar diffusion and reconnection diffusion are very dif- 
ferent. For instance, if the ionization of media is high, am- 
bipolar diffusion predicts that the coupling of mass and mag- 
netic field is nearly perfect. At the same time, reconnection 
diffusion is independent of the ionization but depends on the 
scale of the turbulent eddies and on the turbulent velocities. 
In the paper we explain the physics of reconnection diffu- 
sion both from macroscopic and microscopic points of view, 
i.e. appealing to the reconnection of flux tubes and to the 
diffusion of magnetic field lines. We quantify the reconnec- 
tion diffusion rate both for weak and strong MHD turbulence 
and address the problem of reconnection diffusion acting to- 
gether with ambipolar diffusion. In addition, we provide a 
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criterion for correctly representing the magnetic diffusivity 
in simulations of star formation. We discuss the intimate re- 
lation between the processes of reconnection diffusion, field 
wandering and turbulent mixing of a magnetized media and 
show that the role of the plasma effects is limited to "break- 
ing up lines" on small scales and does not affect the rate of re- 
connection diffusion. We address the existing observational 
results and demonstrate how reconnection diffusion can ex- 
plain the puzzles presented by observations, in particular, the 
observed higher magnetization of cloud cores in comparison 
with the magnetization of envelopes. We also outline a pos- 
sible set of observational tests of the reconnection diffusion 
concept and discuss how the application of the new concept 
changes our understanding of star formation and its numeri- 
cal modeling. Finally, we outline the differences of the pro- 
cess of reconnection diffusion and the process of accumu- 
lation of matter along magnetic field lines that is frequently 
invoked to explain the results of numerical simulations. 



1 Introduction 

Magnetic flux freezing is a key textbook concept with a huge 
impact on astrophysical theory. The concept was first pro- 
posed by Alfven (1942) whose principle of frozen-in field 
lines has provided a powerful heuristic influencing our under- 
standing of many astrophysical processes such as star forma- 
tion, stellar collapse, evolution of accretion disks, magnetic 
dynamo etc. This principle, however, is not universal, as we 
discuss below, its violation entails important consequences 
for star formation. 

The phenomenon of fast magnetic reconnection is an ex- 
ample of failure of the magnetic freezing concept. Fast mag- 
netic reconnection is a type of reconnection which does not 
depend on resistivity and thus should proceed in a media of 
negligible resistivity. Solar flares provide an example of the 
rapid energy release which would not be possible if mag- 
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netic fields were perfectly frozen in (see Yamada, Kulsrud & 
Ji 2010 and references therein). 

The issue of "flux freezing" violation by magnetic recon- 
nection has been known for some time, but was not taken 
very seriously due to the unclear nature of fast reconnection 
(see Zweibel & Yamada 2009 and references therein). In- 
deed, for years it was considered that fast reconnection re- 
quired some special physical conditions and therefore the 
"flux freezing " is fulfilled everywhere apart from some spe- 
cial zones. 

The present day star formation paradigm has been de- 
veloped that flux freezing concept holds. In magnetically- 
mediated star formation theory which was founded by the 
pioneering studies by L. Mestel and L. Spitzer (see Mes- 
tel & Spitzer 1956, Mestel 1966) and brought to the level 
of sophistication by other researchers (see Shu, Adams & 
Lizano 1987, Mouschovias 1991, Nakano et al. 2002, Shu 
et al. 2004, Mouschovias et al. 2006). According to it, 
magnetic fields slow down and even prevent star formation 
if the media is sufficiently magnetized. In the assumption of 
flux freezing of the magnetic field in the ionized component, 
the change of the flux to mass ratio happens due to neutrals 
which do not feel magnetic field directly, but only through 
ion-neutral interactions. In the presence of gravity, neutrals 
get concentrated towards the center of the gravitational po- 
tential while magnetic fields resist compression and there- 
fore leave the forming protostar (e.g. Mestel 1965). This 
makes star formation inefficient for magnetically dominated 
(i.e. subcritical) clouds. The low efficiency of star forma- 
tion corresponds to observations (e.g. Zuckerman & Evans 
1974), which is usually interpreted as a strong argument in 
support of the above scenario. This does not solve all the 
problems as, at the same time, for clouds dominated by grav- 
ity, i.e. supercritical clouds, this scenario does not work as 
magnetic fields do not have time to leave the cloud through 
ambipolar diffusion. Therefore for supercritical clouds mag- 
netic field should be dragged into the star, forming stars with 
magnetizations far in excess of the observed ones (see Galli 
et al. 2006, Johns-Krull 2007). 

What we described above is the initial stage of star for- 
mation. However, magnetic fields are important for other 
stages of star formation, e.g. formation of the accretion 
disks around forming stars. This, as we discuss below, is 
also problematic if one relies on ambipolar diffusion. In 
fact, classical ideas of star formation based exclusively on 
ambipolar diffusion have been challenged by observations 
(Troland & Heiles 1986, Shu et al. 2006, Crutcher et al. 
2009, 2010, see Crutcher 2012 for a review). While the in- 
terpretation of particular observations is the subject of scien- 
tific debates (see Moschovias & Tassis 2009), it is sugges- 
tive that there may be additional processes that the classical 
theory does not take into account. The primary suspect is 
turbulence, which is ubiquitous in the interstellar media and 
molecular clouds (see Larson 1981, Armstrong et al. 1994, 
Elmegreen & Falgarone 1996, Lazarian & Pogosyan 2000, 



Stanimirovic & Lazarian 2001, Heyer & Brunt 2004, Padoan 
et al. 2006, 2009, Chepurnov & Lazarian 2010, Burkhart 
et al. 2010). Turbulence has revolutionized the field of star 
formation (see Vazquez-Semadeni et al. 1995, Ballesteros- 
Paredes et al. 1999, Elmegreen 2000, 2002, McKee & Tan 
2003, Elmegreen & Scalo 2004, MacLow & Klessen 2004, 
McKee & Ostriker 2007) but the treatment of the turbulent 
magnetic fields stayed within the flux freezing paradigm. 

The understanding of flux freezing in turbulent astrophysi- 
cal environments has been challenged relatively recently and 
not all the consequences of this radical change have been 
evaluated so far. Lazarian & Vishniac (1999, henceforth 
LV99) identified magnetic field wandering, which is inher- 
ent property of magnetized turbulent plasma, as the cause of 
fast, i.e. independent of resistivity, magnetic reconnection. 
They showed that in turbulent fluids magnetic fields should 
undergo constant reconnection and change their identity all 
the time. This implies that magnetic fields are not any more 
frozen into a perfectly conducting fluid if this fluid is tur- 
bulent as was explicitly stated first in Vishniac & Lazarain 
(1999). Later, the challenge to the concept of "flux freezing" 
came from another side, i.e. from more formal mathematical 
studies of magnetic fields properties in turbulent fluids (see 
Eyink201 la). Eyink, Lazarian & Vishniac (201 1, henceforth 
ELV1 1) showed the consistency of these two approaches and 
established the equivalence of the LV99 treatment with that 
in more recent mathematical papers. 

While the idea that turbulence can change the reconnec- 
tion rates has been discussed in a number of earlier papers, 
the LV99 model was radically different from its predecessors. 
For instance, Mathaeus & Lamkin (1985, 1986) performed 
2D numerical simulations of turbulence and provided argu- 
ments in favor of magnetic reconnection getting fast. How- 
ever, the physics of the processes that they considered was 
very different from that in LV99 (see more §5.3). 

The LV99 analytical predictions have been successfully 
tested in Kowal et al. (2009) which made it important to 
study astrophysical applications of the model. Magnetic re- 
connection was treated in LV99 for both collisional and col- 
lisionless turbulent plasmas and was extended to the partially 
ionized gas in Lazarian, Vishniac & Cho (2004). This moti- 
vated Lazarian (2005) to identify fast reconnection of mag- 
netic field as an essential process of removing magnetic flux 
at different stages of star formation. Later, in Lazarian et al. 
(2010) this process was termed "reconnection diffusion". 

One may claim that the reconnection diffusion concept 
extends the concept of hydrodynamic turbulent diffusion to 
magnetized fluids. The physics of it is very different from 
the "magnetic turbulent diffusivity" idea discussed within the 
theories of kinematic dynamo (see Parker 1979). Reconnec- 
tion diffusion, unlike "magnetic turbulent diffusivity", deals 
with dynamically important magnetic fields, e.g. with sub- 
Alfvenic turbulence. Thus magnetic fields are not passively 
mixed and magnetic reconnection plays a vital role for the 
process. Recently, the consequences of reconnection diffu- 
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sion have been studied numerically for the diffuse interstellar 
media, molecular clouds and accretion disks (Santos-Lima 
et al. 2010, 2012). These studies should be understood in 
the appropriate context. For example, it is wrong to view 
numerical calculations in Santos-Lima et al. (2010, 2012) 
as the actual justification of the reconnection diffusion con- 
cept. The concept of reconnection diffusion can only be jus- 
tified appealing to the LV99 model and to high resolution 
numerical testing of the latter. Indeed, the testing of the 
reconnection rates predicted in LV99 is performed at much 
higher resolution (see Kowal et al. 2009, 2012, Lazarian 
et al. 2011) than the resolution of the cores and accretion 
disks in our simulations in Santos-Lima et al. (2010, 2012). 
In other words, the numerical effects are under control in 
the simulations focused to test LV99 theoretical predictions, 
while the LV99 theory justifies why one should not be too 
much worried about small scale numerical effects present in 
the astrophysically-movitated set ups exploring astrophysical 
consequences of reconnection diffusion. 

The major shortcoming of our previous papers on recon- 
nection diffusion is that the physical picture of the process 
has never been clearly explained. This causes confusion and 
prevents the heuristic use of the reconnection diffusion con- 
cept. The situation is aggravated by the fact that the textbook 
picture of reconnection involves magnetic fluxes of opposite 
polarities getting into contact. This is clearly different from 
what one expects to see during the star formation process 
where magnetic fields without large scale reversals are be- 
ing dragged together towards the center of the gravitational 
potential. 

The goal of this paper is to clarify what are the actual 
foundations of the reconnection diffusion concept, present 
the physical picture of the diffusion of magnetic field and 
plasmas in the presence of turbulence, to provide the com- 
parison of the predictions based on the reconnection diffu- 
sion concept and observations as well as to present testable 
predictions. We attempt to provide foundations for the al- 
ternative picture of the star formation process based on the 
process of reconnection diffusion. 

In what follows, we discuss the problem of magnetic flux 
removal for star formation in §2, present the description of 
MHD turbulence in §3, describe the model of magnetic re- 
connection in turbulent media in §4, introduce the concept 
of reconnection diffusion in §5. In §6 we discuss the micro- 
scopic picture of reconnection diffusion and in §7 consider 
a possibility of two distinct regimes of reconnection diffu- 
sion. §8 outlines the limitations of numerical studies and 
presents the numerical results of reconnection diffusion for 
different stages of star formation; in §9 we discuss how the 
existing observational data corresponds to the reconnection 
diffusion model and provide observational predictions. Ad- 
ditional consequences of reconnection diffusion for observa- 
tions and numerics are outlined in §10, while the discussion 
and the conclusions are presented in §11 and §12, respec- 
tively. 



2 Star Formation and Magnetic Flux Problem 

2.1 Turbulence in magnetized interstellar plasmas 

A paradigm shift, a concept popularized by the science histo- 
rian Thomas Kuhn (1962), usually delivers a major advance 
in our perception of reality. In their 2004 ARA&A review, 
Scalo & Elmegreen write: "One of the most important devel- 
opments in the field of interstellar gas dynamics during the 
last half century was renewed perception that most processes 
and structures are strongly affected by turbulence. This is a 
paradigm shift unparalleled in many other fields of astron- 
omy, comparable perhaps to the discovery of extrasolar plan- 
ets and cosmological structure at high redshift." 

It has been known for decades that interstellar medium 
(ISM) is driven by violent supernovae explosions (McKee & 
Ostriker 1977). By now it has been accepted that the ISM is 
turbulent on scales ranging from AUs to kpc (see Armstrong 
et al. 1995, Elmegreen & Scalo 2004, Lazarian 2009). Fig. 1 
shows the turbulent power density plotted against the inverse 
of the scale length, with data at large scales, i.e. at small 
wavenumbers q expanded using the Wisconsin H Q Mapper 
(WHAM) data on electron density fluctuations (Chepurnov 
& Lazarian 2010). 

A more direct evidence comes from the observations of 
spectral lines. Apart from showing non-thermal Doppler 
broadening (see Larson 1981), they also reveal spectra of su- 
personic turbulent velocity fluctuations when analyzed with 
techniques like Velocity Channel Analysis (VCA) of Veloc- 
ity Coordinate Spectrum (VCS) developed (see Lazarian & 
Pogosyan 2000, 2004, 2006, 2008) and applied to the obser- 
vational data (see Padoan et al. 2004, 2009, Chepurnov et al. 
2010). 

We should clarify that turbulence in astrophysical envi- 
ronments is the most natural and expected phenomenon. 
Magnetized astrophysical plasmas generally have very large 
Reynolds numbers due to the large length scales involved, as 
well as the fact that the motions of charged particles in the 
direction perpendicular to magnetic fields are constrained. 
Laminar plasma flows at these high Reynolds numbers R = 
VLf/v, where V and Lf are the velocity and the scale of 
the flow, v is fluid viscosity, are prey to numerous linear and 
finite-amplitude instabilities, from which turbulent motions 
readily develop. The drivers of turbulence include super- 
novae explosions that shape the interstellar medium (McKee 
& Ostriker 1977, Nakamura et al. 2006), accretion flows (see 
Klessen & Hennebelle 2010), magneto-rotational instability 
in the galactic disk (Sellwood & Balbus 1999), thermal in- 
stability (see Kritsuk & Norman 2002, Koyama & Inutsuka 
2002), collimated outflows (see Nakamura & Li 2007) etc. In 
fact, present day understanding of interstellar medium views 
molecular clouds as part of interstellar turbulent cascade, i.e. 
as density fluctuations in a compressible flow (see Ostriker, 
Stone, & Gammie 2001, McKee & Ostriker 2007 and ref- 
erences therein). As molecular clouds are magnetized, one 
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Fig. 1. Turbulence in the interstellar gas as revealed by electron 
density fluctuations. "Big Power Law in the Sky" in Armstrong 
et al. (1995) extended using WHAM data. The slope corresponds 
to that of Kolmogorov turbulence. Modified from Chepurnov & 
Lazarian (2010). 



must understand the diffusion of magnetic fields in turbulent 
flows in order to describe star formation. 

2.2 Star formation and role of turbulence 

Star formation is known to be a rather inefficient process. In- 
deed, the mass of the Milky way is Mmw ~ 10 9 solar mass. 
For the typical density of the gas of 50 cm~ 3 the free fall 
time is r ff (3tt jilGpfl 2 « 6 x 10 6 years which provides 
a "natural" star formation rate Mmw / T ff of 200 solar mass 
per year. At the same time the measured star formation rate 
is w 1.3 solar mass per year (Murray & Rahman 2010). It 
is because of this inefficiency we still have large amounts of 
present in interstellar media of spiral galaxies. A traditional 
way of explaining this inefficiency is to appeal to magnetic 
forces preventing gravitational collapse (see Mouschovias & 
Spitzer 1976). 

The role of magnetic fields has been subjected to scrutiny 
from the very beginning of research on star formation the- 
ory (see Mestel & Spitzer 1956). Eventually, star formation 
happens due to gravity, but if the portion of gas is strongly 
magnetized, the magnetic field may prevent such a collapse. 



It is possible to show that if the ratio of the magnetic flux to 
mass is larger than the critical one: 



($/Af) t 



1.8 x 10 3 gauss cm 2 g 1 



(1) 



magnetic field prevents the cloud collapse (see Draine 201 1). 
These strongly magnetized clouds are termed subcritical. 
Such clouds cannot collapse unless than lose their magnetic 
flux. The typical values of magnetization of diffuse gas cor- 
responds to subcritical magnetization. Therefore the com- 
pression of matter together with magnetic field would result 
in subcritical clouds. Naturally, if most of the clouds are sub- 
critical, star formation is inefficient. 

To address the problem of the magnetic field diffusion 
both in the partially ionized ISM and in molecular clouds, 
researchers usually appeal to the ambipolar diffusion con- 
cept (see Mestel 1956, Shu 1983). The idea of the ambipolar 
diffusion is very simple and may be easily exemplified in the 
case of gas collapsing to form a protostar. As the magnetic 
field is acting on charged particles only, it does not directly 
affect neutrals. Neutrals move under the gravitational pull 
but are scattered by collisions with ions and charged dust 
grains that are coupled with the magnetic field. The result- 
ing flow dominated by the neutrals does not drag the mag- 
netic field lines and those will diffuse away dragging through 
the infalling matter. This process of ambipolar diffusion be- 
comes faster as the ionization ratio decreases and therefore, 
becomes more important in poorly ionized cloud cores. The 
corresponding theory of star formation based on ambipolar 
diffusion is well developed (see Mouschovias et al. 2006) 
and it predicts low star formation rates in agreement with ob- 
servations. 

If the cloud is not magnetically dominated, i.e. it is super- 
critical, the above discussion is not applicable and the grav- 
ity should induce a collapse dragging magnetic flux with en- 
trained matter into the forming star on the timescales less 
than the ambipolar diffusion time. This presents twofold 
problem. First of all, simple estimates show that if all the 
magnetic flux is brought together with the material that col- 
lapses to form a star in molecular clouds, then the mag- 
netic field in a proto-star should be several orders of magni- 
tude higher than the one observed in stars (this is the "mag- 
netic flux problem", see Galli et al. (2006) and references 
therein). For instance, T-Tauri stars have magnetic field 
s» 2 x 10 3 Gauss (see Johns-Krull 2007), which amounts to 
($/A/) s» 3 x 10~ 8 gauss cm 2 g~\ which is a million times 
smaller that the flux to mass ratio estimated for one solar 
mass clump in a cloud of density 10 4 cm -3 (see Draine 
201 1). To avoid the "magnetic flux problem" one has to iden- 
tify ways of efficient magnetic flux removal. 

As we discuss earlier, turbulence is an essential process 
for interstellar media. Traditionally, in the textbooks the role 
of turbulence was mostly limited to affecting the virial mass. 
Numerical simulations have shown that turbulence can play 
the dominant role for the formation of the molecular clouds 
(Ballestros-Paredes et al. 2007). Moreover, simulations were 
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indicating the ability of turbulence to change the flux to mass 
ratio even without ambipolar diffusion. Those were inter- 
preted as suggestive of a scenario in which compressible 
turbulence collects matter along magnetic field lines and in- 
duces supercritical star formation (Vazquez-Semadeni et al. 
2011). One might claim that this approach does not re- 
quire magnetic flux diffusion provided that the collapse is 
strictly one-dimensional. As we discuss in §11.6 this sce- 
nario is very restrictive and in this paper we propose a dif- 
ferent solution of the magnetic flux problem for star forma- 
tion, namely, we claim that the turbulent scenario should be 
extended by allowing the process of reconnection diffusion 
that is inevitably induced by MHD turbulence. 

There have been attempts to "turbocharge" ambipolar dif- 
fusion by combining its action with turbulence. We briefly 
discuss these ideas in § 1 1 .4. Our claim in the paper that in 
the presence of turbulence the removal of magnetic flux is 
fast and ambipolar diffusion is not required. 

2.3 Problems of the ambipolar diffusion paradigm 

The alternatives to the classical star formation theory have 
emerged in the last decade as numerical simulations showed 
that star formation can also be slow due to feedback intro- 
duced by turbulence (see McKee & Ostriker 2007, Pudritz 
& Kevlahan 2012). This feedback can disperse clouds be- 
fore they can give birth to new stars. The fact that clouds 
and cores must be necessarily subcritical has also been chal- 
lenged both numerically (see Padoan et al. 2004) and ob- 
servationally (see Crutcher 2012 and references therein). In 
fact, in his 2012 ARA&A review Crutcher states: "There 
is no definitive evidence for subcritical molecular clouds or 
for ambipolar diffusion driven star formation." Below we list 
some other examples where researchers claim that ambipolar 
diffusion is not adequate to explain observations. 

Shu et al. (2006) explored the accretion phase in low- 
mass star formation and concluded that ambipolar diffusion 
could work only under rather special circumstances like, for 
instance, considering particular dust grain sizes. Instead they 
proposed a solution which is based on the magnetic flux dif- 
fusion by Ohmic resistivity. However, to do the job, they 
postulated that the Ohmic resistivity should be increased by 
about 4 orders of magnitude. We feel that there is no physical 
justification for such a resistivity enhancement (see more in 
§11.5). 

Observations of different regions of the diffuse ISM com- 
piled by Troland & Heiles (1986) indicate that magnetic 
fields and density are not correlated. Such a correlation 
would be expected in a naive picture of turbulence where 
compressions of magnetic fields are accompanied by com- 
pressions of density. Magnetic diffusion can explain the ab- 
sence of correlations, but the ambiplar diffusion in diffuse 
ISM is negligible due to high degree of gas ionization. In 
addition, observations of magnetized cores by Crutcher et 
al. (2010) contradict to the predictions of the ambipolar dif- 



fusion paradigm. All these cases are troublesome from the 
point of ambipolar diffusion paradigm but, as we discuss fur- 
ther in the paper (see §9), are consistent with the reconnec- 
tion diffusion concept that we advocate. 

3 Turbulent astrophysical media and its description 

To quantify reconnection diffusion the quantitative descrip- 
tion of astrophysical turbulence is required. 

3.1 Magnetized turbulence in astrophysical plasmas 

In addition to being turbulent, astrophysical plasmas are 
magnetized (see Spitzer 1978, Draine 2011). The magneti- 
zation of astrophysical fluids most frequently arises from the 
dynamo action to which turbulence is an essential component 
(see Schekochihin et al. 2007). In fact, it has been shown that 
turbulence converts in weakly magnetized conducting fluid 
from five to ten percent of the energy of the cascade into 
the magnetic field energy (see Cho et al. 2009). This fraction 
does not depend on the original magnetizationj and therefore 
magnetic fields will come to equipartitiorH with the turbulent 
motions within a few eddy turnover times. 

We deal with magnetohydrodynamic (MHD) turbulence 
which provides a correct fluid-type description of plasma tur- 
bulence at large scales (see §4.2). Astrophysical turbulence 
is a direct consequence of large scale fluid motions expe- 
riencing low friction. The Reynolds numbers are typically 
very large in astrophysical flows as the scales are large. As 
magnetic fields decrease the viscosity for the plasma motion 
perpendicular to their direction, Re numbers get really astro- 
nomically large. For instance, Re numbers of 10 10 are very 
common for astrophysical flow. For so large Re the inner 
degrees of fluid motion get excited and a complex pattern of 
motion develops. 

The drivers of turbulence, e.g. supernovae explosions in 
the interstellar medium, inject energy at large scales and then 
the energy cascades down to small scales through the hier- 
archy of eddies spanning up over the entire inertial range. 
The famous Kolmogorov picture (Kolmogorov 1941) corre- 
sponds to hydrodynamic turbulence, but, as we discuss fur- 
ther, a qualitatively similar turbulence also develops in mag- 
netized fluids/plasmas. Therefore both turbulence and mag- 
netic fields should be dealt with while addressing the prob- 
lem of star formation. Direct observations in Milky Way and 
nearby galaxies provide ample evidence that the star forma- 
tion happens in magnetized turbulent clouds. Due to the pro- 
cess of turbulent dynamo magnetic energy increases fast and 

'This makes the problem of the initial or seed magnetic field, 
that for a long time has worried researchers, rather trivial. Very 
weak magnetic fields, e.g. generated by Bierman battery (see Lazar- 
ian 1992) can be amplified fast in a turbulent plasmas. 

2 In supersonic flows compressibility effects induce deviations 
from the equipartition. 
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one may argue that the formation of stars at high redshifts, 
including the first starts in early universe (see Norman, Wil- 
son & Barton 1980, Abel et al. 2002, Nakamura & Umemura 
2001, Scheicher et al. 2010) could also take place in turbu- 
lent magnetized environments. 

3.2 Strong and weak Alfvenic turbulence 

For the purposes of reconnection diffusion that we de- 
scribe below, Alfvenic perturbations are vital. Numerical 
studies in Cho & Lazarian (2002, 2003) showed that the 
Alfvenic turbulence develops an independent cascade which 
is marginally affected by the fluid compressibility. This ob- 
servation corresponds to theoretical expectations of the Gol- 
dreich & Sridhar (1995, henceforth GS95) theory that we 
briefly describe below (see also Lithwick & Goldreich 2001). 
In this respect we note that the MHD approximation is widely 
used to describe the actual magnetized plasma turbulence 
over scales that are much larger than both the mean free path 
of the particles and their Larmor radius (see Kulsrud 1983, 
2005 and references therein). More generally, the most im- 
portant incompressible Alfenic part of the plasma motions 
can described by MHD even below the mean free path but on 
the scales larger than the Larmor radius (see also §4.2). 

While having a long history of ideas, the theory of MHD 
turbulence has become testable recently due to the advent of 
numerical simulations (see Biskamp 2003) which confirmed 
(see Cho & Lazarian 2005 and references therein) the pre- 
diction of magnetized Alfvenic eddies being elongated in the 
direction of magnetic field (see Shebalin, Matthaeus & Mont- 
gomery 1983, Higdon 1984) and provided results consistent 
with the quantitative relations for the degree of eddy elonga- 
tion obtained in GS95. 

The hydrodynamic counterpart of the MHD turbulence 
theory is the famous Kolmogorov (1941) theory of turbu- 
lence. In the latter theory energy is injected at large scales, 
creating large eddies which correspond to large Re num- 
bers and therefore do not dissipate energy through viscos- 
it)0 but transfer energy to smaller eddies. The process con- 
tinues untill the cascade reaches the eddies that are small 
enough to dissipate energy over eddy turnover time. In the 
absence of compressibility the hydrodynamic cascade of en- 
ergy is ~ vf /r casCt i = const, where vi is the velocity at the 
scale I and the cascading time for the eddies of size I is 
i~cask,i ~ l/vi. From this the well known relation vi ~ Z 1 / 3 
follows. 

A frequent mental picture that astrophysicists have of the 
Alfvenic turbulence is based of Alfven waves with wavevec- 
tors along the magnetic field. This is not true for the 

3 Reynolds number Re = LfV/v = (V / L f ) / \v / L 2 f ) which is 
the ratio of an eddy turnover rate T~ ddy — V/Lf and the viscous dis- 
sipation rate t^ s = r)/ 'L 2 j. Therefore large Re correspond to neg- 
ligible viscous dissipation of large eddies over the cascading time 
Tease which is equal to T e ddy in Kolmogorov turbulence. 



strong Alfvenic turbulence which, similar to its hydrody- 
namic counterpart, can be described in terms of eddieO 
However, contrary to Kolmogorov turbulence, in the pres- 
ence of dynamically important magnetic field eddies become 
anisotropic. At the same time, one can imagine eddies mix- 
ing magnetic field lines perpendicular to the direction of 
magnetic field. For the latter eddies the original Kolmogorov 
treatment is applicable resulting in perpendicular motions 

1/3 

scaling as vi ~ lj_ , where l± denotes eddy scales measured 
perpendicular to magnetic field. These mixing motions in- 
duce Alfvenic perturbations that determine the parallel size 
of the magnetized eddy. The key stone of the GS95 the- 
ory is critical balance, i.e. the equality of the eddy turnover 
time l±/vi and the period of the corresponding Alfven wave 
~ 1\\/Va, where L is the parallel eddy scale and Va is the 
Alfven velocity. Making use of the earlier expression for v\ 

2/3 

one can easily obtain hi ~ , which reflects the tendency 
of eddies to become more and more elongated as the energy 
cascades to smaller scales (see Beresnyak, Lazarian & Cho 
2005). 

It is important to stress that the scales l± and L are mea- 
sured in respect to the system of reference related to the di- 
rection of the local magnetic field "seen" by the eddy. This 
notion was not present in the original formulation of the 
GS95 theory and was added to it in LV99. The local system 
of reference was later used in numerical studies in Cho & 
Vishniac (2000), Maron & Goldreich (2001), Cho, Lazarian 
& Vishniac (2002) testing GS95 theory. In terms of mixing 
motions, it is rather obvious that the free Kolmogorov-type 
mixing is possible only in respect to the local magnetic field 
of the eddy rather than the mean magnetic field of the flow. 

While the arguments above are far from being rigorous 
(see more in a review by Cho & Lazarian 2005), they cor- 
rectly reproduce the basic scalings of magnetized turbulence 
when the velocity equal to Va at the injection scale L. The 
most serious argument against the picture is the ability of ed- 
dies to perform mixing motions perpendicular to magnetic 
field. Jumping ahead of our story, we can mention that this 
ability is related to the ability of magnetic field lines to recon- 
nect fast, i.e. at the rate independent of the fluid resistivity 
(see more in §5.1). 

GS95 theory assumes the isotropic injection of energy at 
scale L and the injection velocity equal to the Alfven ve- 
locity in the fluid Va, i.e. the Alfven Mach number Ma = 
(Vl/Va) = 1, where Vl is the injection velocity. Thus it 
provides the description of transAlfvenic turbulence. This 
model was later generalized for both subAlfvenic, i.e. Ma < 
1, and superAlfvenic, i.e. Ma > 1, cases (see Lazarian & 
Vishniac 1999 and Lazarian 2006, respectively; see also Ta- 
ble 1). Indeed, if Ma > 1, instead of the driving scale L for 

4 The description in terms of interacting wavepackets or modes 
is also possible with the corresponding wavevectors tending to get 
more and more perpendicular to the magnetic field as the cascade 
develops. 
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Table 1 



Regimes and ranges of MHD turbulence 



Type 


Injection 


Range 


Motion 


Ways 


of MHD turbulence 


velocity 


of scales 


type 


of study 


Weak 


V L <V A 


\L , /trans] 


wave-like 


analytical 


Strong 










subAlfvenic 


V l <Va 


[Itrans 5 Imin ] 


eddy-like 


numerical 


Strong 










superAlfvenic 


V l >Va 




eddy-like 


numerical 



L and l m in are injection and dissipation scales 

hrans and I a are given by Eq. J24b and Eq. l !17b . respectively. 



one can use another scale, namely Ia (see Eq. (TPTl l. which 
is the scale at which the turbulent velocity equals to Va- For 
Ma ^> 1 magnetic fields are not dynamically important at 
the largest scales and the turbulence at those scales follows 
the isotropic Kolmogorov cascade vi ~ I 1 / 3 over the range of 
scales At the same time, if Ma < 1, the turbulence 

obeys GS95 scaling (also called "strong" MHD turbulence) 
not from the scale L, but from a smaller scale hrans given 
by Eq. d24l i. while in the range [L,l tr ans] the turbulence is 
"weak". We discuss more superAlfvenic and subAlfvenic 
turbulence in §5.1. 

The properties of weak and strong turbulence are rather 
different. Weak turbulence is wave-like turbulence with wave 
packets undergoing many collisions before transferring en- 
ergy to small scales. It corresponds well to the mental picture 
of turbulence of weakly interacting Alfvenic waves that used 
to dominate astrophysics textbooks. Weak turbulence, unlike 
the strong one, allows an exact analytical treatment (Gaultier 
et al. 2000). On the contrary, the strong turbulence is eddy- 
like with cascading happening similar to Kolmogorov turbu- 
lence within roughly an eddy turnover time. The strong in- 
teractions between wave packets prevent the use of perturba- 
tive approach and do not allow exact derivations. It were the 
numerical experiments that proved the predicted scalings for 
incompressible MHD turbulence (see Cho & Vishniac 2000, 
Maron & Goldreich 2001, Cho et al. 2002, Beresnyak & 
Lazarian 2010, Beresnyak 2011) and for the Alfvenic com- 
ponent of the compressible MHD turbulence|f| (Cho & Lazar- 
ian 2002, 2003, Kowal & Lazarian 2010). 

One also should keep in mind that the notion "strong" 
should not be associated with the amplitude of turbulent mo- 
tions but only with the strength of the non-linear interaction. 
As the weak turbulence evolves, the interactions of wave 
packets get stronger making the turbulence strong. In this 
case, the amplitude of the perturbations can be very small. 



5 For compressible MHD turbulence simulations in Beresnyak et 
al. (2005) and Kowal, Lazarian & Beresnyak (2007) demonstrated 
that the density spectrum becomes more shallow and isotropic as 
the Mach number increases. 



While there are ongoing debates whether the original 
GS95 theory must be modified to better describe MHD turbu- 
lence, we believe that, first of all, we do not have compelling 
evidence that GS95 is not adequate^. Moreover, the proposed 
additions to the GS95 model do not change the nature of the 
physical processes that we discuss below. 



4 Magnetic reconnection of turbulent magnetic field 

In what follows we describe the model of reconnection of 
magnetic field proposed in LV99. 

4.1 Turbulence in reconnection zone 

Magnetic reconnection is a fundamental process that vio- 
lates the frozen-in state of magnetic flux. We would like to 
stress that we are discussing the case of dynamically impor- 
tant magnetic field, including the case of a weakly turbulent 
magnetic field. The case of a weak magnetic field which can 
be easily stretched and bent by turbulence at any scale up to 
the dissipation one is rather trivial and of little astrophysical 
significance for star formation (see more in §5.1). Indeed, 
at sufficiently small scales magnetic fields get dynamically 
important even for superAlfvenic turbulence (see Eq. (fTTTi). 

Within the picture of eddies mixing perpendicular to the 
local magnetic field (see §3.2), it is suggestive that magne- 
tized eddies can induce turbulent diffusion similar to ordi- 
nary hydrodynamic eddies. This is a rather counter-intuitive 
notion in view of the well-entrenched idea of flux being 
frozen in to astrophysical fluids. As it is explained in ELV1 1 
the frozen-in condition is not a good approximation for tur- 



Recent work by Beresnyak & Lazarian (2009a, 2010) shows 
that present day numerical simulations are unable to reveal the ac- 
tual inertial range of MHD turbulence making the discussions of 
the discrepancies of the numerically measured spectrum and the 
GS95 predictions rather premature. In addition, new higher reso- 
lution simulations by Beresnyak (2011) reveal the predicted —5/3 
spectral slope. 
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Sweet-Parker model 




Being more quantitative, one can write 
A 

Vrec =VAJ—, 
J-'x 



(2) 




blow up 



Fig. 2. Upper panel: Sweet-Parker reconnection. A is limited by re- 
sistivity and small. Middle panel: reconnection according to LV99 
model. A is determined by turbulent field wandering and can be 
large. Lower panel: magnetic field reconnect over small scales. 
From Lazarian, Vishniac & Cho (2004). 



meaning that v rec <C Va if A <C L x . There are different ways 
to derive Sweet-Parker formulae for the reconnection rate 
(see Parker 1979). One way is to consider the Ohmic diffu- 
sion of magnetic field lines (see ELV11). The mean-square 
vertical distance that a magnetic field-line can diffuse by re- 
sistivity in time t is 

(y 2 (t))~\t, (3) 

where A = ijc 2 /Air is Ohmic diffusivity. The field lines are 
advected out of the sides of the reconnection layer of length 
L x at a velocity of order Va- Thus, the time that the lines 
can spend in the resistive layer is the Alfven crossing time 
tA = Lx/Va- Thus, field lines can only be merged that are 
separated by a distance 

A = y/(y 2 (t A )) ~ v / A*7 = L X /VS, (4) 



bulent fluidfl The violation of the frozen in condition also 
follows from LV99 model of reconnection (see discussion in 
Vishniac & Lazarian 1999). 

A picture of two flux tubes of different directions which 
get into contact in 3D space is a generic framework to de- 
scribe magnetic reconnection. The upper panel of Figure 
|2] illustrates why reconnection is so slow in the textbook 
Sweet-Parker model. Indeed, the model considers magnetic 
fields that are laminar and therefore the frozen-in condition 
for magnetic field is violated only over a thin layer domi- 
nated by plasma resistivity. The scales over which the re- 
sistive diffusion is important are microscopic and therefore 
the layer is very thin, i.e. A <C L x , where L x is the scale at 
which magnetic flux tubes come into contact. The latter is 
of the order of the diameter of the flux tubes and typically 
very large for astrophysical conditions. During the process 
of magnetic reconnection all the plasma and the shared mag- 
netic flu)H arriving over an astrophysical scale L x should be 
ejected through a microscopic slot of thickness A. As the 
ejection velocity of magnetized plasmas is limited by Alfven 
velocity Va, this automatically means that the velocity in the 
vertical direction, which is reconnection velocity, is much 
smaller than Va- 



'Formal mathematical arguments on how and why the frozen-in 
condition fails may be found in Eyink (201 1) (see also Eyink 2007, 
2009). 

8 Figure[2]presents a cross section of the 3D reconnection layer. 
A shared component of magnetic field is present in the generic 3D 
configurations of reconnecting magnetic flux tubes. 



where S is Lundquist number, 

S = L X V A /X- (5) 

Combining Eqs. (0 and © one gets the famous Sweet- 
Parker reconnection rate, 

V rec ,SP = VA/VS. (6) 

The LV99 model generalizes the Sweet-Parker one by ac- 
counting for the existence of magnetic field line stochastic - 
ity (Figure [2] [lower panels]). The shown turbulence is sub- 
Alfvenic and the mean field is clearly defined. At the same 
time turbulence induces magnetic field wandering. This wan- 
dering was quantified in LV99 and it depends on the intensity 
of turbulence. The vertical extend of wandering of magnetic 
field lines that at any point get into contact with the field of 
the other flux tube was identified in LV99 with the width of 
the outflow region. In other words, the LV99 model of recon- 
nection makes use of the fact that in the presence of magnetic 
field wandering, which is a characteristic feature of magne- 
tized turbulence in 3D, the outflow is no more constrained by 
the narrow resistive layer, but happens through a much wider 
area A defined by wandering magnetic field lines. An im- 
portant consequence of this is that as turbulence amplitude 
increases, the outflow region and therefore reconnection rate 
also increases, which entails the ability of reconnection to 
change its rate depending on the level of turbulence. The 
latter is important both for understanding the dynamics of 
magnetic field in turbulent flow and for explaining flaring re- 
connection events, e.g. solar flares. 
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We should note that the magnetic field wandering is mostly 
due to Alfvenic turbulence^ and the corresponding expres- 
sions for A arising from the field wandering were obtained 
in LV99. An alternative derivation of the A was obtained 
in analogy with the Sweet-Parker derivation above in ELV1 1 
appealing to the concept of Richardson (1926) diffusion. 

Richardson diffusion (see Kupiainen et al. 2003) implies 
the mean squared separation of particles 



(\x 1 (t)-x 2 (t)\'' 



(7) 



where t is time, e is the energy cascading rate and () de- 
note an ensemble averaging. For sub Alfvenic turbulence 
e « u\/ {VaL), where ml is the injection velocity and L is 
an injection scale (see LV99) and therefore analogously to 
Eq. (|4]l one can write 



et 3 A ^L x (L x /L)^ 2 M_ 



(8) 



where it is assumed that L x < L. Combining Eqs. (O and 
(O one recovers the LV99 expression for the rate of magnetic 
reconnection (see also ELV1 1) 



iLV99 ^V A (L x /L) 1 / 2 M_ 



2 

A- 



(9) 



in the limit of L x < L. Analogous considerations allow to 
recover the LV99 expression for L x > L, which differs from 
Eq. (O by the change of the power 1/2 to— 1/2. 

It is important to stress that Richardson diffusion ulti- 
mately leads to the diffusion over the entire width of large 
scale eddies once the plasma has moved the length of one 
such eddy. The precise scaling exponents for the turbulent 
cascade does not affect this result, and all of the alternative 
scalings considered in LV99 yield the same behavior. 

The predictions of the turbulent reconnection rates in 
LV99 were successfully tested 3D numerical simualtions in 
Kowal et al. (2009) (see also Lazarian et al. 2010 for an 
example of higher resolution runs). In Figure [3] we see the 
results for varying amounts of input power, for fixed resistiv- 
ity and injection scale as well as for the case of no turbulence 
at all. The line drawn through the simulation points is for 
the predicted scaling with the square root of the input power. 
The agreement between equation and Figure[3]is encour- 
aging. Similarly the dependences of the reconnection rate on 
the injection scale and on the Ohmic and anomalous resistiv- 
ity were successfully tested. 

The testing of LV99 predictions provides more confidence 
to the theory and stimulates to think of its applications, e.g. 



As discussed in LV99 and in more details in ELV11 the mag- 
netic field wandering, turbulence and magnetic reconnection are 
very tightly related concepts. Without magnetic reconnection, prop- 
erties of magnetic turbulence and magnetic field wandering would 
be very different. For instance, in the absence of fast reconnection, 
the formation of magnetic knots arising if magnetic fields were not 
able to reconnect would destroy the self-similar cascade of Alfvenic 
turbulence. The rates predicted by LV99 are the rates required to 
make Goldreich-Sridhar model of turbulence self-consistent. 
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Fig. 3. Reconnection speed versus input power for the driven tur- 
bulence. We show the reconnection speed plotted against the input 
power for an injection wavenumber equal to 8 (i.e. a wavelength 
equal to one eighth of the box size) and a resistivity v u . The dashed 
line is a fit to the predicted dependence of P 1 / 2 (see eq. (3)). The 
horizontal line shows the laminar reconnection rates for each of the 
simulations before the turbulent forcing started. Here the uncer- 
tainty in the time averages are indicated by the size of the symbols 
and the variances are shown by the error bars. Modified from Kowal 
et al. (2009). 



for star formation. One should keep in mind that the LV99 
model assumes that the magnetic field flux tubes can come at 
arbitrary angle, which corresponds to the existence of shared 
or guide field within the reconnection layeiB 

4.2 LV99 Reconnection and Plasma effects 

For years plasma effects have been considered essential for 
fast magnetic reconnection (see Shay et al. 1998, Daughton 
et al. 2006, 2008). On the contrary, LV99 makes use of 
the MHD approximation. To issue of the justification of the 
MHD description of plasmas in the LV99 model was recently 
revisited in ELV11. One can think of three relevant length- 
scales: the ion gyroradius pi, the ion mean-free-path length 
£mfp,i> an d the scale L s of large-scale variation of magnetic 
and velocity fields. Astrophysical plasmas can be "strongly 
collisional" if £ m fp,i *C Pi, and can be described as fluids. 
The interiors of stars and accretion disks present examples of 
such plasmas. Another case is "weakly collisional" t m jp,i 3> 
pi plasmas. The ratio of the mean free path to the gyroscale 

A v A 
InA c ' 



Pi 



(10) 



follows from the expression for the Coulomb collision fre- 
quency (see Fitzpatrick 201 1, eq.(1.25)), where A = AttuX'jj 



The model in LV99 is three dimensional, and it is not clear 
to what extend it can be applied to 2D turbulence (see discussion 
in ELV11 and references therein). However, the cases of pure 2D 
reconnection and 2D turbulence are of little practical importance. 
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is the plasma parameter, or number of particles within the 
Debye screening sphere. Hot and rarified astrophysical plas- 
mas are weakly coupled which entails A being large, for in- 
stance, of the order of 10 9 or more for the warm compo- 
nent of the interstellar medium or solar wind (see Table 1 
in ELV11). For such ratio the expansion over small ion gy- 
roradius results in "kinetic MHD equations". Those differ 
from the standard MHD by having anisotropic pressure ten- 
sor (see more discussion in Kowal et al. 201 la and references 
therein). 

In addition, plasmas that are not strongly collisional can 
be divided into two subclasses: "collisionless" plasmas for 
which the mean free pass is larger than the largest scales 
of interest £ m fp.i ^ L s „ and "weakly collisional" plasmas 
for which the opposite is true, i.e. L s ^> l m fp,i- In the lat- 
ter case the"kinetic MHD" description allows further reduc- 
tion in complexity at scales larger the mean free path l m f p ^. 
This, as discussed in ELV11, reproduces a fully hydrody- 
namic MHD description at those scales. For instance, the 
warm ionized ISM is "weakly collisional", while the solar 
wind interaction with the magnetosphere is "collisionless." 

Additional simplifications are possible when (a) turbulent 
fluctuations are small compared to the mean magnetic field, 
(b) have length-scales parallel to the mean field much larger 
than perpendicular length-scales, and have frequencies low 
compared to the ion cyclotron frequency. This set of as- 
sumptions (a), (b) and (c) is adopted in the GS95 theory 
and for "gyrokinetic approximation" (see Schekochihin et 
al 2007). At length-scales, i.e. at scales much larger than 
pi, an important simplification takes place. At those scales 
incompressible shear- Alfven wave modes have get indepen- 
dent from the compressive modes and described by the sim- 
ple "reduced MHD" (RMHD) equations (see GS95, Cho & 
Lazarian 2003). This fact is of major importance for the 
LV99 justifying the use of the analysis based on an incom- 
pressible MHD fluid model. 

For most astrophysical problems, including the problem of 
star formation, one has deal with the reconnection at length 
scales much larger than ion gyroradius, i.e. L s ^> pi. For 
the solar wind at in the Earth magnetosphere where some of 
the reconnection measurement have been done (Retino et al. 
2007), the ion gyroradius pi w 6 x 10 6 cm, i.e. comparable 
with L s . In these circumstances, for the span of scales from 
pi to the electron gyroradius p e the plasmas is described by 
an ion kinetic equation and electron "reduced MHD" (ER- 
MHD) equations for kinetic Alfven modes (see Schekochi- 
hin et al. 2009). This is the domain of Hall physics, with ion 
and electrons moving very differently and "Hall reconnec- 
tion" being expected (see Uzdensky & Kulsrud 2006). Mag- 
netic field wandering still can take place due to the electron 
fluid (see Daughton et al. 2011) but the nature of the tur- 
bulent cascade is different (see Cho & Lazarian 2004, 2009, 
Schekochihin et al. 2007, 2009). We hope that future re- 
search will clarify to what extend features of the LV99 model 
carries over to this very different regime of reconnection in 



3D: ?:■>{>' ■'. ;;..■> k. -A, c> r K ' 




n nn r - 

16 18 20 22 24 26 28 30 

Time 




Fig. 4. Upper p/o/:Reconnection rate in the presence of anomalous 
resistivity and turbulence. Changes of the anomalous resistivity do 
not change the reconnection rate. Lower plot: Visualization of re- 
connection region. Turbulence is weakly driven and the big changes 
in magnetic field lines are due to reconnection. Plotted with data in 
Kowal et al. (2009). 



magnetosphere. 

Within the LV99 model, the reconnection rate is deter- 
mined by large scale magnetic wandering (see Figure|2]for a 
pictorial representation and FigureQT]for the analysis of sim- 
ulations), while small scale plasma effects are irrelevant for 
the global reconnectior0 at least in the fully ionized plasma 
(see Lazarian et al. 2004). This conclusion is supported by 
simulations in Kowal et al. (2009) where plasma effect were 
simulated by using anomalous resistivity, i.e. the resistiv- 
ity that depends on the value of the current. Figure [4] shows 
that substantial variations of the anomalous resistivity do not 
change the resulting reconnection rates. Note that the numer- 
ical effects produce also a sort of anomalous resistivity on 

"it is shown in LV99 that if the reconnection is calculated as- 
suming that the bottle neck is due to Ohmic resistivity, the recon- 
nection rate gets much larger than the Alfven speed. This is due to 
the effect of 3D turbulent magnetic fluxes getting in contact over 
many independent patches. The total cumulative rate of reconnec- 
tion "cutting magnetic field lines" becomes very large and it is the 
outflow of magnetized plasma from the reconnection region that 
enters to limit the overall reconnection speed. Naturally, increas- 
ing the local reconnection speed of magnetic patches either due to 
plasma effects or, alternatively, numerical effects in computer sim- 
ulations does not increase the reconnection speed of the large scale 
turbulent magnetic fluxes. 
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the scales comparable to the grid size. Therefore our results 
testify that, in the presence of turbulence, numerical effects 
do not dominate simulations in terms of reconnection. We 
should mention that although plasma effects do not change 
the global reconnection rate, they can be important for other 
processes, e.g. the acceleration of electrons (see §10.1). 

4.3 Alternative views on fast reconnection 

Alternative models of magnetic reconnection appeal to dif- 
ferent physics to overcome the limitations of the Sweet- 
Parker model. In the Petschek (1964) model of reconnection 
the reconnection layer opens up to enable the outflow which 
thickness does not depend on resistivity. There the extend 
of the current sheet gets of the order of microscopic A and 
therefore Eq. (0 provides v rec w Va- For years this had been 
considered the only way of make magnetic reconnection fast. 
To realize this idea inhomogeneous resistivity, e.g. anoma- 
lous resisitivity associated with plasma effects, is required 
(see Shay & Drake 1998, Shay et al. 1998, 1999, Bhattachar- 
jee et al. 2005, Cassak et al. 2006, 2008, 2009). However, 
for turbulent plasmas, the effects arising from modifying the 
local reconnection events by introducing anomalous resistiv- 
ity are negligible as confirmed e.g. in Kowal et al. (2009). 
Other effects, e.g. formation and ejection of plasmoids (see 
Shibata & Tanuma 2001, Lourreiro et al. 2007, Uzdensky et 
al. 2010, Huang et al. 2011) may be important for initially 
laminar environments and contribute to the onset of recon- 
nection. However, for the problems of star formation, it is 
more important to account for the pre-existing turbulence for 
which LV99 model is directly applicable. 

Turbulence is known to make accelerated diffusive pro- 
cesses. Therefore it is not surprising that it had been ap- 
pealed as the way of speeding up the reconnection prior to 
LV99 study. Nevertheless, LV99 model is radically differ- 
ent from its predecessors which also appealed to the effects 
of turbulence. For instance, unlike Speiser (1970) and Ja- 
cobson (1984) the model does not appeal to changes of the 
microscopic properties of the plasma. The closest in its ap- 
proach to LV99 among papers dealing with effects of turbu- 
lence was the work of Matthaeus & Lamkin (1985, 1986) 
who studied the problem numerically in 2D MHD and who 
suggested that magnetic reconnection may be fast due to a 
number of turbulence effects, e.g. multiple X points and tur- 
bulent EMF However, the physics of reconnection discussed 
in Matthaeus & Lamkin (1985, 1986) is very different, as 
they did not realize the key role played by magnetic field- 
line wandering, which is the corner stone idea of the LV99 
model. They did not obtain a quantitative prediction for the 
reconnection rate, as did LV99. 

A number of earlier papers may be seen as indirect evi- 
dence of fast reconnection in turbulent fluids. For instance, a 
study of tearing instability of current sheets in the presence 
of background 2D turbulence and the formation of large- 
scale, long-lived magnetic islands has been performed in 



dPolitano et al.Lll989t) . They present evidence for "fast en- 



ergy dissipation" in 2D MHD turbulence and show that this 
result does not chan ge as they change the reso lution. A more 
recent study of of dMininni & Pouquet 1 120091) also provides 
evidence for "fast dissipation" but in 3D MHD turbulence. 
This phenomenon is consistent with the idea of fast recon- 
nection, but cannot be treated as a direct evidence of the pro- 
cess. It is very clear that fast energy dissipation and fast 
magnetic reconnection are distinclty different physical pro- 
cesses. In addition, a paper by Galsgaard and Nordlund, 
(IGalsgaard & Nordlund L 1 1997a) . might also be interpreted 



as providing indirect evidence for fast reconnection. The au- 
thors noted that in their simulations they could not produce 
highly twisted magnetic fields. One of the interpretations of 
this finding could be the relaxation of magnetic field via re- 
connection. In this case, this observations could be related to 
the numerical finding of (?) which shows that reconnecting 
magnetic configurations spontaneously get chaotic and dis- 
sipate, which in its turn may be related to the predictions in 
LV99 (see more in Lapenta & Lazarian 2011). However, in 
view of the uncertainties of the numerical studies, it is diffi- 
cult to be confident of this connection. 

LV99 model deals with balanced turbulence where the en- 
ergy flows in opposite directions are equal. If the latter is 
not true, MHD turbulence is imbalanced, or has non-zero 
cross helicity. Solar wind presents an example of system with 
imbalanced turbulence. There have been recent attempts to 
study reconnection in systems with a flow and imbalanced 
turbulence (Yokoi & Hoshino 2011). We feel that to obtain 
quantitative predictions one needs to use the scaling proper- 
ties of the turbulence as it is done in LV99. At the same time 
the theory of strong imbalanced MHD turbulence is still con- 
troversial. Am ong the existing theories of imbalanced tur- 
bulence (see jL ithwick et al., 2007; Beresnvak & Lazarian , 
l2008t IChandranl l2008UPerez & Boldvrevl 120091). all, but 
(Beresnva k & Lazarianl 2008b seem to contradict to numer- 
ical testing in Beresnyak & Lazarian (2009b). We defer the 
discussion of reconnection in imbalanced turbulence to fu- 
ture publication^ 

It is important to stress that while possible reconnection 
schemes may be numerous, the generic astrophysical recon- 
nection model should satisfy several constraints. A funda- 
mental consideration for such a model is that they must ex- 
plain fast reconnection in collisional and collisionless plas- 



~We expect the effect of field wandering to play the crucial role 
for the reconnection with non-zero cross-helicity. This wandering 
should be determined using the theory of strong imbalanced MHD 
turbulence, for instance, by one in Beresnyak & Lazarian (2008) 
if higher resolution and longer averaging future testings confirm it. 
We also expect that similarly as a successful model of imbalanced 
MHD turbulence must produce GS95 scaling for the case of zero 
cross helicity, the reconnection in flows with imbalanced turbulence 
should converge to LV99 predictions for the zero cross helicity. For 
very high cross helicity field wandering may be small and other, e.g. 
plasma effects, become important. 
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mas. At the same time, to explain flares, it should be possi- 
ble for reconnection to be delayed for significant amounts of 
time. The reconnection model should be able to operate in 
the turbulent environment as astrophysical media are turbu- 
lent. As far as we know, only LV99 model satisfies to all of 
these requirements. 

For instance, LV99 model explains the accumulation of 
flux in highly magnetized plasmas if the level of turbulence is 
low. It is possible that tearing mechanisms may provide orig- 
inal perturbations stimulating the development of turbulence 
and reconnection that it induces. As the outflow within the 
reconnection region gets turbulent, it induces more of field 
wandering of the reconnecting fluxes and therefore higher 
reconnection speed (see more in Lazarian & Vishniac 2009). 
This introduces the positive feedback which results in a flare. 
LV99 predicted that turbulence from neighboring regions can 
also ignite reconnection and the observation of the initiation 
of flares by incoming Alfvenic waves was reported by Sych 
et al. (2009). 

We would like to stress, that magnetic turbulence and re- 
connection are intrinsically connected. Therefore there is no 
"magic" reconnection rate, e.g. O.IVa, that satisfies all the 
requirement^ For instance, to avoid formation of magnetic 
knots turbulent eddies should be able to reform their mag- 
netic field structure over the time scales of their turnover. 
For transAlfvenic turbulence this means the reconnection 
rate ~ Va- This stringent constraint is satisfied by the LV99 
model (see §5.1). 



5 Reconnection diffusion concept 

5.1 Diffusion in magnetized turbulent fluid 

The exact treatment of diffusion in turbulent fluid is rather 
complicated and it requires dealing with magnetic reconnec- 
tion as a part of the process. Therefore we first discuss an 
illustrative example with pure hydro turbulence. 

Hydrodynamic Diffusion 
To illustrate analytical approaches to diffusion in a turbu- 
lent fluid one can consider first incompressible unmagnetized 
fluid where the velocity U is decomposed into a regular part 
V and a fluctuation v. Averaging the Navier-Stockes equa- 
tions one gets 
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where P is the average pressure, (..) denote averaging proce- 
dure and the indices indicate vectors and the standard sum- 
mation convention. The term (viVj) is not specified by the 
equation and its approximation involves the different "clo- 
sures" (see Monin & Yaglom 1975 for a discussion of vari- 



The latter was a sort of Holy Grail for many researchers study- 
ing reconnection. 



ous hydrodynamic closures). A similar problem emerges in 
the description of the diffusion induced by turbulence. 

Consider as an example the diffusion of a scalar quantity 
s, e.g. a passive impurity, e.g. heavy elements in the ISM. If 
only molecular diffusion were present, the rate of transport 
of s in a;-direction would be given by Fick's law, namely, 



ox 



(12) 



where D is the molecular diffusion coefficient. This law 
is being modified by the advection of the quantity s by the 
transport induced by the velocity field. Then one can get the 
advective - diffusion equation (see Fisher 1979): 
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Decomposing the field s into the mean s m and fluctuating 
s / parts one gets after averaging from Eq. Q~3] 
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where the last term corresponds to the transport of the field 
by turbulent fluctuations compared to the Fick's law based 
on molecular diffusion (see Eq. [L21 . The frequently used 
"closure" appeals to "eddy diffusivity" coefficients defined 
as follows 



( V 3 S f) 



ds 
dxj 



(15) 



where, in general, coefficients e may depend on the direc- 
tion. For hydrodynamic turbulence ei ~ vil, which results in 
Richardson (1926) diffusion coefficient I 4 / 3 if velocity 
field is Kolmogorov, i.e. Vi ~ I 1 / 3 . For the diffusion at the 
largest scales induced by turbulence with the injection scale 
L and velocity Vl 



Khydro ~ LVl 



(16) 



corresponds to the maximal diffusivity of unmagnetized tur- 
bulent fluid. The above treatment can approximate the dy- 
namics of fluid perpendicular to the local magnetic field, pro- 
vided that the reconnection is fast. This is the issue that we 
address below. 

Diffusion in Magnetized Fluid 
The introduction of magnetic field complicates the process 
of diffusion as it forces one to account for the back reac- 
tion of magnetic field. The latter was a hotly debated subject 
in the dynamo theory (see Parker 1979). For infinitesimally 
weak magnetic field which backreaction is negligible the tur- 
bulent eddies bend magnetic field lines over all scales up to 
the Ohmic dissipation one and the disparity of scales present 
in the problem of reconnection described in §4.1 does not 
emerge. This way of reasoning resulted in the concept of 
magnetic turbulent diffusivity within kinematic dynamo ac- 
cording to which the diffusion of magnetic field should is 
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similar to the diffusion of passive scalar in hydrodynamic tur- 
bulence. However, this regime of dynamically unimportant 
magnetic field presents a highly unrealistic case of marginal 
astrophysical importance^- 

The problem of diffusion of magnetized fluid when mag- 
netic fields are dynamically important demands that the issue 
of magnetic reconnection is to be addressed. Otherwise, the 
existence of eddies is questionable. Mixing motions in MHD 
turbulence require that reconnection events in MHD turbu- 
lence should happen through every eddy turnover. This is, 
however, what the LV99 model predicts. Indeed, for small 
scales magnetic field lines are nearly parallel and, when 
they intersect, the pressure gradient is not V\/l\\ but rather 
(l±/l?i)Vjl, since only the energy of the component of the 
magnetic field that is not shared is available to drive the out- 
flow. On the other hand, the characteristic length contraction 
of a given field line due to reconnection between adjacent ed- 
dies is This gives an effective ejection rate of V A /l\\- 
Since the width of the diffusion layer over the length h is 
l± Eq.© should be replaced by Vr ~ Va(1±/1\\), which pro- 
vides the reconnection rate Va/Zii, which is just the nonlinear 
cascade rate on the scale In. This ensures self-consistency 
of critical balance for strong Alfvenic turbulence in highly 
conducting fluids (LV99). In fact, if not for the LV99 recon- 
nection the buildup of unresolved magnetic knots would be 
unavoidable, flattening the turbulence spectrum compared to 
the theoretical predictions. The latter contradicts both Solar 
wind measurements and numerical calculation^ 

Let us first consider the maximal rate allowed by the re- 
connection diffusion, i.e. evaluate the diffusivity arising 
from the largest eddies. Dealing with reconnection diffu- 
sion one should consider all regimes of MHD turbulence, su- 
perAlfvenic, transAlfvenic and subAlfvenic. We start with a 
superAlfvenic regime, i.e. Ma > 1. Magnetic field gets dy- 
namically important as soon as its energy density exceeds the 
energy of eddies at the Ohmic dissipation scale, which trans- 
lates into Alfvenic velocity getting larger than the velocity 
of eddies at the Ohmic dissipation scale or the ion Larmor 
radius, whichever is larger. As the velocity in Kolmogorov 
turbulence scale as vi ~ I 1 / 3 , it is clear that even weak mag- 
netic field can make a significant impact on the dynamics of 
the smallest eddies. In view of that it is advantageous to in- 
troduce a scale at which the magnetic field gets dynamically 
important and the nature of the turbulence changes from hy- 



14 In the case of a dynamically unimportant field the magnetic dis- 
sipation and reconnection happens on the scales of the Ohmic diffu- 
sion scale and the effects of magnetic field on the turbulent cascade 
are negligible. However, turbulent motions transfer an appreciable 
portion of the cascading energy into magnetic energy (see Cho et al. 
2009; also §3.1). As a result, the state of intensive turbulence with 
negligible magnetic field is short-lived. 

15 If magnetic reconnection is slow then, as was claimed by Don 
Cox (private communication), the interstellar medium should be- 
have not like a fluid, but more like felt or Jello. 



drodynamic to MHD (see Lazarian 2006), namely, 

l A = L{V A /V L f = LMl* (17) 

where, similar to the notation in the hydrodynamic case, L 
and Vl are the injection scale and injection turbulent veloc- 
ity, respectively. If the mean free path of particles is larger 
than Ia, the scale Ia may act as an effective mean free path 
in terms of particle diffusion along magnetic fields. This is 
an important consideration for the diffusion of heat in colli- 
sionless fluid, but it is not so for the diffusion of magnetic 
field and plasmas induced by turbulence. The corresponding 
diffusion coefficient for the maximal rate coincides with its 
hydrodynamic counterpart given by Eq. ( fT6] i, i.e. 



^supA — ^hydro- 



(18) 



This is due to the fact that for the largest eddies of su- 
perAlfvenic turbulence are marginally affected by magnetic 
field. 

For low Alfvenic Mach numbers, i.e. for V A 3> Vl at large 
scales ~ L the turbulence is weak (see Ng & Bhattachargee 
1997, LV99, Gaultier et al. 2000) and magnetic fields are 
slightly perturbed by propagating Alfven waves. The wave 
packets in weak turbulence evolve changing their perpendic- 
ular scale while their scale l\\ along the magnetic field 
does not change. The diffusion that is being induced by weak 
turbulence is substantially reduced compared to the case of 
hydrodynamic turbulence. It can be estimated as 



(19) 



where d is the random walk of the field line over the wave 
period ~ The weak turbulence at scale L evolves over 
the non-linear evolution time (see review by Cho, Lazarian 
& Vishniac 2003 and references therein) 



T-iVA/VL) 2 ^ 1 



(20) 



The transverse displacement of magnetic field lines in the 
perpendicular direction is a result of random walk 



(y 2 ) ~ {™)d 2 



(21) 



According to LV99, the transverse displacement of magnetic 
field lines over distance x is a spatial random walk given by 
the equation d(y 2 ) jdx ~ L{ Vl / Va) 4 , which results in 



{y 2 )~Lx(V L /V A y 



(22) 



Combining Eqs. <|T9), ©, ED and (ED one gets the diffu- 
sion coefficient for the weak turbulence 



K we ak ~ LVl(V L /V A ) 3 = LV L M$ 



(23) 



which is smaller than its hydrodynamic counterpart by the 
factor M\ < 1. 

The additional contribution to diffusivity in the case 
V A ^> Vl comes from scales at which magnetic turbulence 



14 



A. Lazarian: Reconnection and Star Formation 



gets strong. As scaling of weak turbulence predicts Vi ~ 
Vl(/_l/Z) 1/2 (LV99), at the scale 

kr m .~L(V L /V A ) 2 = LM\ (24) 

the critical balance condition l\\/V A ~ l±/Vl is getting satis- 
fied making turbulence strong. It is easy to see that the ve- 
locity corresponding to l trans is V trans ~ V L (V L /V A ). For 
strong turbulence the diffusion 

l^stronq ~V trans ltrans~LV L {V L /V A r (25) 

which coincides with Eq. d23T i. indicating that the en- 
hanced diffusivity of smaller eddies in strong MHD turbu- 
lence regime can produce as efficient magnetic field mixing 
as the diffusivity induced by weak turbulence at the injection 
scale. This coincidence illustrate the deep connection of the 
weak and strong turbulence in terms of the mixing that these 
processes induce. 

Dealing with subAlfvenic turbulence one should distin- 
guish between the diffusivity parallel and perpendicular to 
magnetic field. Magnetic field in a turbulent fluid changes 
the diffusion of plasma with particles moving along magnetic 
field lines. At the same time, turbulent eddies in the direc- 
tion perpendicular to the local direction of magnetic field, 
according LV99, are similar to the Kolmogorov picture. In- 
deed, in the presence of fast reconnection magnetic mixing 
is not inhibited for motions perpendicular to the local direc- 
tion of magnetic field. As the local direction of the magnetic 
field varies in respect to the mean magnetic field, the diffu- 
sivity induced by mildly subAlfvenic and trans-Alfvenic tur- 
bulence is not expected to be very different for the directions 
parallel and perpendicular to the mean magnetic field. At the 
same time, superAlfvenic turbulence, naturally, induces an 
isotropic diffusion. 

The transAlfvenic case of M A ~ 1 is the case which the 
GS95 model in its original formulation deals with. It should 
not be special and should correspond to substituting M A = 
1 in the expressions obtained for superAlfvenic and sub- 
Alfvenic cases, e.g. in Eq. (l25T l. 

As we mentioned earlier, the diffusion of magnetic field 
and plasmas at scales smaller than the scales of injection L 
as well as the scale transfer to the Alfvenic turbulence ltrans 
for the case of subAlfvenic turbulence, is subject to Richard- 
son diffusion (see Eq. (|7]l). This is the type of diffusion that 
gets accelerated with time as well as with scale involved, i.e. 
"superdiffusion". We discuss below that dealing with star 
formation one is interested in the rate of the diffusion from 
a given clump or a cloud. In this case one can use the scale 
and velocity dispersion of the cloud as the proxies of L and 
Vl- In the case of subAlfvenic turbulence, the diffusion co- 
efficient can be approximated as 

Kcloud ~ V c l ou dlcloud(v cloud /V A ) 3 (26) 

while for the case of turbulence being transAlfvenic or su- 
perAlfvenic the same estimate but without M\ factor should 
be used. 



Before reconnection YZ view After reconnection 




Fig. 5. Motion of matter in the process of reconnection diffu- 
sion. 3D magnetic flux tubes get into contact and after reconnection 
plasma streams along magnetic field lines. Left panel: XY projec- 
tion before reconnection, upper panel shows that the flux tubes are 
at angle in X-Z plane. Right panel: after reconnection. 

5.2 Physical picture of reconnection diffusion in the ab- 
sence of gravity 

The description above provides the mathematical framework 
of turbulent diffusivity in a homogeneous magnetized fluid. 
To what extent these results carry over to the mixing of highly 
inhomogeneous magnetic fields important for star formation 
is being discussed below. 

We should mention that the concept of reconnection dif- 
fusion is based on LV99 model and was first discussed in 
Lazarian (2005) in the context of star formation. The con- 
cept, however, is also applicable to heat transfer (see Lazar- 
ian 2006) as well as the turbulent transfer of heavy elements, 
i.e. "metals", in the galaxy (see Lazarian 201 1). The dynamo 
implications of the LV99 reconnection model were discussed 
in LV99 and we may claim that the reconnection diffusion 
of dynamically important magnetic fields should be used in- 
stead of the "turbulent diffusion" that was introduced to de- 
scribe kinematic dynamo, i.e. the dynamo applicable to the 
case of extremely weak magnetic fields onhl 16 ! (see §5.1). In 
terms of star forming clouds, the most important type of dy- 
namo is the turbulent dynamo (see Cho et al. 2010 and ref. 
therein). Such a dynamo may be essential for the clouds in 
the early Universe when the plasmas were only weakly mag- 
netized. Turbulent dynamo can also play role in increasing 



It is important to stress that fast reconnection does not change 
the helicity of the magnetic flux and therefore it does not solve the 
problem of helicity in astrophysical dynamos (see Vishniac & Cho 
2001). Therefore the simulations where researchers use enhanced 
many orders of magnitude resistivity in an attempt to mimic effects 
of turbulence smoothing on magnetic field are in error. However, 
the transport of magnetic flux and smoothing that does not change 
helicity are important ingredients of any mean field dynamo that 
reconnection diffusion takes care of. 
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magnetization of superAlfvenic turbulent clouds. However, 
within the present paper we are mostly concerned with the 
issues of magnetic field removal via reconnection diffusion. 

If magnetic field lines preserve their identify the diffusion 
of charged particles perpendicular to magnetic field lines is 
very restricted and the mass loading of magnetic field lines 
does not change. However, LV99 model (see also ELV11) 
suggests that the standard assumptions are violated if mag- 
netized fluids are turbulent. As a result, in the presence of 
MHD turbulence, relative motions of plasma perpendicular 
to magnetic field not only possible, but inevitable. 

We shall first illustrate the reconnection diffusion process 
showing how it allows plasma to move perpendicular to the 
mean inhomogeneous magnetic field (see Figure |5). The set 
up is relevant to what we encounter in star formation. Two 
magnetic flux tubes with entrained plasmas intersect each 
other at an angle and due to reconnection the identity of 
magnetic field lines change. Before the reconnection plasma 
pressure P p i as ma m the tubes is different, but the total pres- 
sure Ppiasma + Pmagn is the same for two tubes. This is 
a textbook situation of a stable equilibrium. If plasmas are 
partially ionized then slow diffusion of neutrals which do not 
feel magnetic field directly can make gradually smoothen the 
magnetic field pressure gradients. We claim that in the pres- 
ence of turbulence a different process takes place. 

Magnetic field lines in the presence of turbulence are not 
parallel. Such field lines can reconnect and do reconnect all 
the time in a turbulent flow (see §4. 1). The process of recon- 
nection changes the topology of the initial magnetic configu- 
ration and connects magnetic fields with different mass load- 
ing and plasma pressures. As a result, plasmas stream along 
newly formed magnetic field lines to equalize the pressure 
along flux tubes. Portions of magnetic flux tubes with higher 
magnetic pressure expand as plasma pressure increases due 
to the flow of plasma along magnetic field lines. The entropy 
of the system increases with magnetic and plasma pressures 
becoming equal through the volume. In other words, an effi- 
cient process of the diffusion of plasmas and magnetic field 
takes place and this process does not rely on the partial ion- 
ization of the material. In the absence of gravity, the ef- 
fect of this process is to make magnetic field and plasmas 
more homogeneously distributed. The effective diffusion of 
both magnetic field and plasmas is about LV s t r eam, where 
Vstream arises due to plasma pressure difference along the 
parts of the reconnected flux tube. Therefore it is of the or- 
der of the sound speed. A "microscopic" picture of the same 
process is presented in §6. 

If the process shown in Figure [5] were the only one, the 
speed of reconnection diffusion would be limited by the 
speed of plasma motion along magnetic field lines. This is 
not the case, however. The exchange, as a result of reconnec- 
tion, of parts of flux tubes with entrained material is another 
process essential for the process of reconnection diffusion. 
This ensures, for instance, that for supersonic turbulence, the 
exchange is happening at turbulent velocities. 




Fig. 6. Reconnection diffusion: exchange of flux with entrained 
matter. Illustration of the mixing of matter and magnetic fields due 
to reconnection as two flux tubes of different eddies interact. Only 
one scale of turbulent motions is shown. In real turbulent cascade 
such interactions proceed at every scale of turbulent motions. 



To get a clear mental picture of what is going on, consider 
a toy model of turbulence with only one scale of motions. 
Thus, similar to Figure [5] where instead of considering the 
entire multitude of flux tubes of different sizes we focused on 
what is happening with two flux tubes, here we consider just 
two neighboring eddies as is illustrated by Figure [6] Mag- 
netic flux tubes moving with the eddies reconnect and plas- 
mas and magnetic fields associated with one eddy becomes 
a part of the other. This induces efficient diffusion of both 
magnetic field and plasmas. Indeed, if, as in the case illus- 
trated by Figure [5] the densities of plasma along magnetic 
field lines are different in the two flux tubes, the reconnection 
in Figure [6] creates new flux tube with columns of entrained 
dense and rarefied plasmas. The situation is similar to the 
earlier discussed case with plasma moving along magnetic 
fields and equalizing the pressure within the newly formed 
flux tubes. As a result, eddies with initially different plasma 
pressure exchange matter and equalize plasma pressure. This 
process can be described as the diffusion of plasma perpen- 
dicular to the mean magnetic field. 

In reality, for turbulence with the extended inertial range, 
the shredding of the columns of plasmas with different den- 
sity proceeds at all turbulence scales making speed of plasma 
motion irrelevant for the diffusion. For the case of strong tur- 
bulence the diffusion of matter and magnetic field is given by 
Eq. (23. 

Naturally, the process of reconnection diffusion also takes 
place when the pressure of plasmas is the same throughout 
the volume. This is important for the diffusion of plasma im- 
purities and heat. The mixing is happening as new magnetic 
flux tubes are constantly formed from magnetic flux tubes 
that belong to different eddies. It is clear that plasmas which 
were originally entrained on different flux tubes get into con- 
tact along the flux tubes created through reconnection. The 
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process similar to the depicted one takes place at different 
scales down to the scale of the smallest eddies. At the small- 
est scales, the microscopic diffusivity of plasmas takes over. 

The efficiency of reconnection diffusion depends on the 
scale of the motions involved. The process of reconnection 
diffusion may be illustrated with the case of the diffusion 
of impurity from a blob of the size a (Figure [7J. This set- 
ting allows us to consider homogeneous turbulence, which is 
simplifies the analysis. Turbulence is characterized by its in- 
jection scale L max , its dissipation scale L m i n and its inertial 
range [L m i n ,L max \. Consider Alfvenic eddies perpendicular 
to magnetic field lines. If turbulent eddies are much smaller 
than a, i.e. a 3> L m i n they extend the spot acting in a ran- 
dom walk fashion. For eddies much larg er than the blob, i.e. 
a <C L m i n they mostly advect the spoO If a is the within the 
inertial range of turbulent motions, i.e. L m i n < a < L max 
then a more complex dynamics of turbulent motions is in- 
volved. This is the case of Richardson diffusion (see Eq. 

If the blob a is not just an impurity, but has density dif- 
ferent from the density of the surrounding flow, the process 
of reconnection diffusion depends on the properties of tur- 
bulence within and outside the blob. The level and scale of 
turbulence may differ within and outside the blob. We may 
use the idealized picture in Figure [7J to get a qualitative in- 
sight, nevertheless. If the blob is constrained by gravity, as 
this is relevant to star formation, and is turbulent up to the 
largest scale the diffusivity of magnetic field from the blob 
can be roughly estimated as av a , where v a is the velocity 
at the scale a. This estimate assumes that scale a is within 
the range of strong MHD turbulence. Then the reconnec- 
tion diffusion from the volume scales as a 4 / 3 . If a is within 
the range of weak turbulence the reconnection diffusivity is 
given by Eq. ( 126*1 ). 

Our example in Figure |7J illustrates the diffusion perpen- 
dicular to magnetic field. As we mentioned in §3.2, Alfvenic 
motions that are most efficient in mixing perpendicular to 
the local direction of magnetic field. This direction, in gen- 
eral, does not coincide with the mean magnetic field direc- 
tion. Therefore in the system of reference related to the mean 
magnetic field (and to the external observer) the diffusion of 
magnetic field and plasmas will happen both parallel and per- 
pendicular to the mean magnetic field direction. However, 
the weaker the perturbations of the magnetic field, i.e. the 
smaller the Vl /Va ratio, the more anisotropic the reconnec- 
tion diffusion is. 

5.3 Reconnection diffusion in the presence of gravity 

In the presence of gravitational forces acting upon plasmas, 
the diffusion will also increase the entropy of the system, 
allowing the "heavy fluid", i.e. gas, to be concentrated and 

17 There is also expansion of the spot arising from the Lyapunov 
deviation of the flow lines as we discuss in §6.1. 




Fig. 7. Reconnection diffusion depends on the size a of the zone 
from which the diffusion happens. Different regimes emerge de- 
pending on the relation a to the sizes of maximal and minimal 
eddies present in the turbulence cascade. Eddies perpendicular to 
magnetic field lines correspond to Alfvenic turbulence. The plots 
illustrate heat diffusion for different regimes. Upper plot corre- 
sponds to a being less than the minimal size of turbulent eddies; 
Middle plot corresponds to a being less than the damping scale of 
turbulence; Lower plot corresponds to a within the inertial range of 
turbulent motions. This is the case of Richardson diffusion. 



"light fluid", i.e. magnetic field, to leave the the gravitational 
potential"*!. 

Consider the idealized system of plasmas and magnetic 
field in a uniform directed downwards gravitational field 
with the acceleration g. In the thermodynamic equilib- 
rium the plasmas will have the Boltzmann distribution 
with pexp[— migz/kT e ff], where the effective temperature 
for supersonic turbulence can be roughly estimated from 
kT e ff ~ rridomVl /2, m& om is the mass of dominant species 
in the flow. As a result, the light fluid, namely magnetic 
field, would fill the entire volume and have the same pres- 
sure. Therefore the magnetization of the media measured in 
terms of magnetic flux to plasma mass ratio is the lowest at 
the bottom and highest at the top of the system. 

The redistribution of the magnetic field is induced by the 
reconnection diffusion. However, one should clarify what 
is the diffusion coefficient for the magnetic field that is in- 
duced by the process. We discussed in §5.2 that the diffusion 
happens in a different way at different scales. It is gener- 
ally accepted that the turbulence in molecular clouds is a part 
of a big power law turbulent cascade (see Figur^U and the 
discussion in §3.1). The gravitational cores that get rid of 
the magnetic field excess are much smaller than the scale of 
interstellar turbulent driving, i.e. ~ 100 pc. Assuming that 



The difference of the reconnection diffusion and ambipolar dif- 
fusion is that the former is associated with turbulent motions of mat- 
ter. Therefore, too intensive turbulence may upset the virial balance 
and instead of magnetic field diffusion, a dispersion of the entire 
cloud can take place. 
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clouds are parts of the same turbulent cascade existing in the 
diffuse gas one can use the idealized sketch in Figure [7] to 
estimate the efficiency of reconnection diffusion. It seems 
natural to associate the size of blob a with the virial radius 
of the core r vir , which for the case of the turbulence support 
of the core is GM/r vir ~ V? b , where G is the gravitational 
constant, M is the core mass and V tur b is the turbulent ve- 
locity associated with the cloud at the scale r V i r . 

6 Reconnection diffusion and the identity of magnetic 
field lines 

Below we consider the process of reconnection diffusion mi- 
croscopically, at the level of individual field lines. 

6.1 Explosive diffusion of magnetic field lines in turbu- 
lent flows 

A textbook description of magnetic field lines in perfectly or 
nearly perfectly conducting fluid assumes that the line pre- 
serves its identity. There is a number of ways to see that this 
is impossible in a turbulent fluid (see ELV1 1). 

We shall start by showing that the Richarson diffusion 
(see Eq. (01) produces very non-trivial results. Consider 
the problem of separating particles in Kolmorogov turbu- 
lence. The separation between two particles d/dtl(t) ~ 
v(l) ~ al 1 / 3 , where a is proportional to a cubical root of 
the energy cascading rate. The solution of this equation is 
= \l^ 3 + a(t — to)} - 1 / 3 , which provides Richardson dif- 
fusion^ of I 2 ^t 3 . However, as correctly stressed by Greg 
Eyink (201 1), the odd feature of this solution is that the pro- 
vides this type of fast separation even if the initial separation 
of particles is zero, which means the violation of Laplacian 
determinism. Mathematically the above paradox is resolved 
by accounting to the fact that turbulent field is not differen- 
tiabl^l and therefore the initial value problem does not have 
a unique solution. Thinking physically, we cannot assume 
that the turbulence is present up to Iq = 0. 

Although the previous example dealt with hydrodynamic 
turbulence, the essential features of this example cary over to 
the case of MHD turbulence, as in the plane perpendicular to 
the local direction of magnetic field, strong MHD turbulence 
satisfies the Kolmogorov descriptiorl^l 

"Richardson diffusion presents an example of superdiffusion, i.e. 
diffusion process for which I 2 ^t l3 ,/3> 1. The important conse- 
quences of Richardson diffusion have been studied for heat transfer 
and cosmic ray propagation (see Lazarian 2006, Yan & Lazarian 
2008). 

20 The Kolmogorov velocity field is Holder continuous, i.e. 

|v(ri)-T>(ra)|<C|ri-ra| 1/3 . 

21 This is not only similarity in terms of the spectrum. Cho, Lazar- 
ian & Vishniac (2003) showed that in the local system of reference 
the intermittency of turbulence is also similar to the hydrodynamic 
one. 



The study in LV99 revealed the Richardson-type diffusion 
of magnetic field lines. Figure [8] illustrates the loss of the 
Laplacian determinism for magnetic field lines. In analogy 
with the illustrative example above, the final line spread l± 
does not depend on the initial separation of the field lines. 
This is a remarkable effect that provides a microscopic pic- 
ture of reconnection diffusion based on the description of 
magnetic field lines rather than on the reconnection of well 
organized flux tubes in §5.1. 

We shall consider tracing magnetic field lines in the real- 
istic turbulence with the dissipation scale l m in,X> where, as 
everywhere in the paper, _L denotes the scale perpendicular to 
the local magnetic field. We feel that this avoids some of the 
paradoxes discussed in ELV1 1 and also allows to treat a more 
generic case of astrophysical turbulence. However, we first 
briefly stress a couple of points presented in ELV11. First 
of all, resistivity, whatever its nature, introduces stochastic 
forcing in the description of magnetic field line dynamics. In- 
deed, the induction equation with the resistive term r/AB sig- 
nifies stochasticity associated with Ohmic diffusion. There- 
fore the definition of the magnetic field line on scales affected 
by resistivity is not deterministic. In addition, one has to ac- 
cept that the magnetic field line motion is a concept defined 
by convention and not testable experimentally. This point 
is stressed in the literature (see Newcomb 1958, Vasylianas 
1972, Alfven 1976, ELV11), but sometimes forgotten. Thus 
magnetic field lines may be tagged by ions that start at the 
same field line (see Figure |9|l In the case of smooth lami- 
nar magnetic field and ideal MHD equations the motions of 
ions will reveal magnetic field lines and two ions on the same 
field line will always remain on the same line. The situation 
is radically different in the presence of turbulence. 

In our thought experiment we shall trace ions moving with 
the same velocity and separated perpendicular to magnetic 
field by a distance of a Larmor radius po. As we discuss 
further, if the separation is less rhoo than that one can ap- 
peal to plasma effects to increase the separation to po- Let 
us assume that the minimal scale of turbulence l m in,± > Po- 
In this situation the dynamics of ions can be approximated 
by the dynamics of charged particles in toy model of "a sin- 
gle scale MHD turbulence" discussed in Rechester & Rosen- 
bluth (1978). Indeed, the turbulent motions at the critically 
damped scale l m in,± are dominant for shearing and steer- 
ing matter and magnetic field on the smaller scales. The 
Rechester & Rosenbluth (1978) theory predicts the Lypunov 
growth of the perpendicular separation of ions, i.e. the sep- 
aration gets poexp(l/l m i n 11, where I is the distance traveled 
by ions and l m i n 11 is the scale parallel scale of the critically 
damped eddies with the perpendicular scale l m i n ,± (see also 
Narayan & Medvedev 2001, Lazarian 2006). Thus to get 
separated by the distance l m in,X the ions should travel the 
so-called Rechester-Rosenbluth distance: 

LRR ~ lmin,\\l n (lmin,±. /po) (27) 
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Fig. 8. Particle tracing magnetic field lines may start at different ini- 
tial locations shown as coaxial ellipsoids. However after a period of 
time the field line spread over a larger volume and the final position 
of the field lines does not correlate with their initial position. 

which at most a dozen times larger than the microscopic 
scale l m i n .\\. As soon as ions get separated over the distance 
of l m in,± they get into different eddies and the process of 
Richardson diffusion starts. Thus after a relatively short pe- 
riod when ions move in a correlated manner remembering 
their original position, a stochastic regime when the initial 
vicinity of the ions is completely forgotten takes over. As we 
used the ions as tracers of magnetic field, we can talk about 
the stochasticity of magnetic field lines (as traced by ions). 

If the minimal turbulence scale is equal to p > the argu- 
ments above only get stronger, as from the very beginning 
the ions may experience stochastic turbulence driving and get 
uncorrelated. In many cases, for instance, in fully or mostly 
ionized ISM (e.g. with the ionization larger than 93% as dis- 
cussed e.g. in Lazarian, Vishniac & Cho 2004), the Alfvenic 
turbulent cascade at po gets continued as a whistler cascade 
involving only electrons. Such cascading provides stochas- 
ticity below po. With the whistler scaling as vi ~ Z 2 / 3 (see 
Cho & Lazarian 2009 and references therein) the "whistler- 
induced" Richardson diffusion should go as I 2 ~ t 6 , inducing 
fast separation of magnetic field lines that can be now traced 
by electrons and inducing stochastic perturbations on ion tra- 
jectories. This only makes the case for the stochasticity when 
the ions are separated initially over less than po more evident. 

Consider the case of ions separate by less than the po dis- 
tance and turbulence truncated at scales larger than pq. The 
generalized Ohms law can be written as 

I „ JxB 

E = — u x B + 7/j_J^ +77||J|i H 

, c nec 

VP e m e f dJ ■. , 1 _,A ,„„. 
- + — | ^- + V(uJ+Ju JJ) , (28) 

nec ne z \ot ne J 




Fig. 9. Upper plot Ions tracing the same magnetic field line. The 
diffusion and decorrelation arises from plasma or Ohmic effects as 
well Rechester-Rosenbluth effect. Middle plot Ions separated by 
scales much larger than the ion Larmor radius are further separated 
by the Rechester-Rosenbluth effect. Lower At scales larger than 
the turbulence damping scale the Richardson diffusion takes over 
resulting in explosive separation of field lines. 



when quasi-neutrality and m e <C m., are assumed (see Bhat- 
tacharjee et al. 1999). The electric fields on the righthand 
side are, respectively, the motional field, Ohmic fields as- 
sociated to perpendicular and parallel resistivities, the Hall 
field, a contribution from the electron pressure tensor, and 
electron inertial contributions. All of these terms contribute 
to the slippage of magnetic field-lines through the plasma on 
the scales of the order of po. This introduces small scale dif- 
fusion and therefore stochasticity for the ions initially very 
close in their position (ELV1 1). Therefore the ions are bound 
to separate fast however close they were initially. The diffu- 
sion introduced by the terms in the Ohm's law is important 
on the scales po, but for the reconnection diffusion the exact 
form of these terms does not matter. The initial separation 
of the particles is being fast forgotten and it is the large scale 
turbulence that determines the macroscopic stochasticity of 
magnetic field lines. 

For the purposes of discussing reconnection diffusion, we 
have chosen ions as our trace particles. One, however, may 
argue that electrons as current carrying agents are more ap- 
propriate. This does not change our arguments above, how- 
ever, as using the Larmor radius of electron peiectron instead 
of po in Eq. d27l ) changes the result by an insignificant factor. 
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6.2 Spontaneous stochasticity of magnetic field lines 
and reconnection diffusion 

The process of magnetic field lines becoming stochastic in 
turbulent fluids is also called "spontaneous stochasticity" 
(see Eyink 201 1, ELV1 1). We believe that the way of tracing 
magnetic field lines with ions discussed in §6.1 has a very 
transparent meaning in terms of diffusion of plasmas in as- 
trophysical conditions. The common wisdom underlying the 
star formation research, for instance, is based on the picture 
of laminar field lines where nearby ions stay all the time en- 
trained on the same field line provided that the non-ideal ef- 
fects, e.g. resistivity, are negligibly small. It is evident from 
§6.1 that this is not correct for realistically turbulent fluids. 

To proceed with our discussion of the physics of reconnec- 
tion diffusion, consider two turbulent volumes at a distance 
between them (see Figure [TUb. Each of the volumes has its 
own set of magnetic field lines. However, as the field lines 
wander and spread as the consequence of the Richardson dif- 
fusion they overlap in the volume Ai nt and their identity as 
associated with the particular volume is lost (Figure[lOl left). 
Magnetic field lines reconnect and the newly formed lines al- 
low plasma exchange between the initially disconnected vol- 
umes. If plasma and magnetic field pressures in the volumes 
was different, one can easily seen that this picture based on 
the diffusion of magnetic field lines is analogous to the pic- 
ture describing the exchange of plasmas between magnetic 
flux tubes in Figure [5] 

Figure \W[ left, illustrates the spread of magnetic field 
lines in the perpendicular direction as magnetic field lines 
are traced by particles moving along them. Using Eq. ( l22b 
one can get the RMS separation of the magnetic field lines 
(LV99) 



Sli 



L WaJ 



(29) 



i.e. l\ is proportional to i|. This regime identified in LV99 is 

a Richardson diffusion but in terms of magnetic field linesF^l 
The numerical testing of this prediction of LV99 is shown in 
Figure QT| This is the case of diffusion in space. Eq. d29l ) 
allow one to calculate the distance In at which the magnetic 
field lines of regions separated by l± start overlapping^ i.e. 
l\\ « (Sl' 2 ± L)^ 3 (V A /V L )^ 3 . Thus for sufficiently large l {] all 
parts of the volume of magnetized plasmas get connected. 
In other words, the entire volume becomes accessible to par- 
ticles moving along magnetic field lines. Naturally, in this 
situation the customary for the star formation community no- 
tion of flux to mass ratio loses its original sense. 



~ This regime induces perpendicular "superdiffusion" in terms of 
cosmic rays and other charged particles streaming along magnetic 
field lines. 

23 Incidentally, in terms of reconnection, Eq. (??) expresses the 
thickness of outflow denoted as A in Figure [2] Substituting this 
value in Eq. (0 one recovers the LV99 reconnection rate, i.e. Eq. 
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Fig. 10. Microscopic physical picture of reconnection diffusion. 
Magnetized plasma from two regions is spread by turbulence and 
mixed up over Aj n t. Left panel: Description of the process in terms 
of field wandering in space. Right panel. Description of the spread 
in time. 



In addition, in turbulent plasmas, our turbulent volumes 
spread due to Richardson diffusion with 51^ ~ t 3 . This pro- 
cess is illustrated by Figure \TU\ right. The magnetized plas- 
mas spreads by subAlfvenic turbulence over a larger vol- 
ume, while motions of plasma along magnetic field lines are 
ignored. The two initially disconnected volumes get over- 
lapped and the magnetic field and plasma of two regions gets 
mixed up. Again, one easily can see that the process that we 
described in terms of magnetic field lines is similar to the 
one we described in terms of reconnected magnetic field flux 
tubes in Figure [6] This is expected, as magnetic reconnec- 
tion is an intrinsic part of the picture of MHD turbulence that 
governs the dynamics of magnetic field lines (LV990 

For the partially ionized gas the following questions are 
practically important. Does the Rechester-Rosenbluth length 
given by Eq. ( |27] i present a serious bottleneck for Richar- 
son diffusion on the larger scales? Does the Richarson diffu- 
sion provide faster diffusion rates compared to the ambipolar 
diffusion? The latter question is addressed in § 1 1 .4 where 
we argue that if ambipolar diffusion is too efficient, we ex- 
pect turbulence to decay, while if strong MHD turbulence is 
present in the volume this means that the reconnection diffu- 
sion dominates. The answer to the former question is more 
subtle. There has not yet been a study of the Richardson dif- 
fusion in a partially ionized gas. However, appealing to the 
equivalence of fast reconnection and spontaneous stochastic- 
ity revealed in ELV1 1, we can appeal to the work by Lazar- 
ian et al. (2004) that demonstrates that the reconnection in 
a partially ionized gas is fast. This justifies our application 
of the Richardson diffusivity to the partially ionized gas. In- 
tuitively, one can argue that the Richardson diffusion on the 



One can also claim that spontaneous stochasticity of magnetic 
field in turbulent fluids is an underlying process that governs mag- 
netic reconnection (ELV11). As we mentioned earlier, fast recon- 
nection makes MHD turbulence theory self-consistent. 
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scales larger than the damping scale makes the latter irrel- 
evant similar to the case of damping scale in fully ionized 
plasma that we discussed earliea 25 l 

Naturally, in the presence of diffusion, ions that we use to 
trace magnetic field spread in the volume, reproducing the 
results of spreading of magnetic field in the turbulent vol- 
ume via reconnection diffusion that we argued in §5.2. In 
the presence of gravity, it is obvious that the lighter fluid of 
magnetized ions should tend to escape the gravitational po- 
tential as we discussed in §5.3. In other words, tracing of 
ions provides us with the microscopic picture of reconnec- 
tion diffusion. 

6.3 Magnetic field wandering and reconnection diffu- 
sion 

Magnetic field wandering described in LV99 is a implemen- 
tation of Richardson diffusion (ELV11). This very wander- 
ing has then been then used in the literature to describe the 
properties of cosmic ray diffusion and heat transport (see 
Narayan & Medvedev 2001, Lazarian 2005, Yan & Lazarian 
2008) and has been tested numerically (Maron et al. 2004, 
Lazarian et al. 2004). Figure [TT] shows numerical results on 
magnetic field separation. The regime corresponding to the 
Richardson-type scaling is clearly seen at the scales corre- 
sponding to the inertial range of the turbulence. At the scales 
larger than the injection scale the usual diffusion regime 
takes over. The latter regime was discussed decades ago 
within earlier models of MHD turbulence (see, for instance, 
classical papers by Parker 1965, Jokipii 1973). 

Neither of the regimes of magnetic wandering can be un- 
derstood as static with magnetic field lines preserving their 
identify in turbulent flows. Magnetic field lines should 
constantly reconnect inducing the exchange of plasmas and 
magnetic field inducing the reconnection diffusion process 
that we advocate in this paper. In view of this, interesting 
ideas on the non-trivial nature of magnetic field lines can be 
found in a prophetic book by Parker (1979). Nevertheless, 
Parker (1979) does not formulate the reconnection diffusion 
or spontaneous stochasticity concepts. 

6.4 MHD and plasma-based descriptions of reconnec- 
tion diffusion 

In the previous section we proved that LV99 ensures efficient 
diffusion of matter and magnetic fields in turbulent fluids. 
One may wonder whether other models of fast reconnection, 
e.g. collisionless reconnection, can also induce reconnection 
diffusion. 



25 In fact, Lazarian et al. (2004) show that the turbulence in par- 
tially ionized gas demonstrates several regimes, including the in- 
termittent "resurrection" of turbulence cascade at scales less that 
the ambipolar damping scale. These effects make magnetic fields 
stochastic on scales much less that the naively estimated l m i n ^\\ that 
enters Eq. ((27}. 
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Fig. 11. Magnetic field wandering from numerical simulations. 
Separation of magnetic field lines shown. Modified from Lazarian 
et al. 2004. 



First of all, we would like to stress that the LV99 model 
is not in conflict with the studies of magnetic reconnection 
in collisionless plasmas that have been a major thrust of the 
plasma physics community (see Shay et al. 1998, Daughton 
et al. 2006). Unlike latter studies, LV99 deals with turbulent 
environments. It shows, as we discussed in §4. 1 that local re- 
connection rates are influenced by plasma effects, e.g. kinetic 
effects of Hall effects, but the overall or global reconnection 
rate, i.e. the rate at which magnetic flux tubes reconnect, is 
determined by the turbulent broadening of the reconnection 
region. Thus, in the turbulent astrophysical media the rate of 
reconnection is not going to be affected by additional mech- 
anisms. As a result, the reconnection diffusion will proceed 
with its maximal rate limited by turbulent motions only. 

If, on the contrary, the media that we deal with is not tur- 
bulent, reconnection diffusion does not take place even in the 
presence of fast reconnection. Turbulent mixing is a neces- 
sary condition for reconnection diffusion to exist. If other 
mechanisms of reconnection induce turbulence, this turbu- 
lence will induce reconnection diffusion and we return back 
to the case above. 

Plasma effect might be potentially important on a more 
subtle level, however. As field wandering is essential for re- 
connection diffusion, one may wonder to what extend the 
plasma effects can be neglected while describing field wan- 
dering. Indeed, as we mentioned in §6.1 the non-ideal terms 
in the Ohms law d28l > provide the stochasticity of charge 
carriers and therefore of the magnetic field lines that these 
charges trace. Consider, for instance, Hall term which is 
most commonly invoked in the literature of fast reconnec- 
tion. Hall term was also invoked in a studies of magnetic field 
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loss by circumstellar accretion disk^l (see Krasnopolsky et 
al. 2010). The usual criterion for the Hall term to dominate 
is that the electron flow velocity is dominated by the cur- 
rent. However, to dominate magnetic field stochasticity the 
criterion should be different as the correlations of the Hall 
velocity are short ranged. Assuming the small-scale equipar- 
tition of velocity and magnetic field, B(r) 2 /4ir ~ pv(r) 2 /2 
and turbulent correlation of velocities (viVj) ~ Cr 1 , one gets 
for the Hall velocity Vho.ii ~ J jne — cV x B /Airne correla- 
tions 

(V Ha ii,iVH a u,j) ~ (^^) 2 A(B l B 3 ) (30) 

The right hand sight Eq. d30l can be estimated as 
f , c j AiTnpCr 1 ^ 2 , which is much smaller than the corre- 
lation of the turbulent velocities if the distance between point 
of correlation r 3> c 2 rrii/ Anne 2 = Si, where Si is the ion in- 
ertial skin depth. This estimate is consistent with a more 
elaborate one in ELV1 1. Therefore even large Hall velocities 
do not affect meandering of magnetic field lines on the scale 
much larger that the inertial skin depth. The reconnection 
diffusion applicable to star formation deals with scales 3> Si. 

To finish with the discussion of plasma effects, we should 
mention that the Hall MHD (HMHD) is frequently presented 
as a proper way to describe reconnection in astrophysical 
systems. However, it is shown in ELV11 that HMHD is 
rarely applicable to the actual astrophysical plasmas. In- 
deed, the derivation of Hall MHD based on collisionality re- 
quires that the ion skin-depth Si must satisfy the conditions 
Si 3> S ^> l m fp,i> where S is the scale of large-scale vari- 
ations of magnetic field. The second inequality is needed 
so that a two-fluid description is valid at the scales of in- 
terest, while the first inequality is needed so that the Hall 
term remains significant at those scales. However, substitut- 
ing Si = pi j \fpl into the expression for the Coulomb colli- 
sional frequency (see Eq. ( [Tol l) yields the result 

— T^-oc— . (31) 

Si InA c 

where A = AnpuX^ is the number of particles in the Debye 
sphere. For weakly coupled astrophysical plasmas A is really 
large (see table in EVL1 1) and therefore £ m fp,i 3> Si, unless 
the ion temperature is extremely low. Thus, Hall MHD is 
valid only for cold, dense plasmas,e.g. that produced by the 
MRX reconnection experiment (e.g. Yamada et al. 2010), 
but not in the conditions of the diffuse ISM and molecular 
clouds where star formation takes place. 

7 Hypothetical weak regime of reconnection diffusion 

Our discussion above shows that, similar to the case of hy- 
drodynamic turbulence, strong MHD turbulence induces ef- 

26 To see the effect the authors had to adopt the Hall term much 
larger than its value for the adopted parameters of the media. 



ficient diffusion. In this regimes the initial separation of par- 
ticles and magnetic field lines does not affect their final po- 
sition after a sufficiently extended period of time (see our 
illustration of the Richardson diffusion in §6.1). For super- 
Alfvenic and transAlfvenic turbulence the regime of acceler- 
ated diffusion spans from the scale of dissipation to the one of 
injection. For subAlfvenic turbulence the Richardson regime 
covers scales from the dissipation scale to the scale of the 
transition to the strong turbulence, i.e. t tra ns given by Eq. 
d24i >. The latter is true provided that the dissipation scale is 
larger than ttrans- 



What would happen if the turbulence dissipation scale is 
larger than t tr ans^ This can result in a very different regime 
of diffusion. Indeed, the diffusion induced by weak tur- 
bulence obeys the ordinary diffusion laws with the squared 
displacement proportional to time, i.e. 8 2 ~ t, rather than 
i 3 for the Richardson diffusion. The trajectories of parti- 
cles exponentially depart, but they remember their starting 
point. Therefore the reconnection is expected to depend on 
the plasma microphysics and will not be fast. This is a very 
special regime of very weak driving and very strong dissi- 
pation. This regime can be called "weak reconnection diffu- 
sion". In this situation, one would expect interesting new ef- 
fects, e.g. the joined action of ambipolar diffusion and weak 
reconnection diffusion. 



We feel that the situation with weak reconnection diffusion 
requires further studies and no definitive conclusions are pos- 
sible at the moment. Potentially, weakly driven turbulence in 
a partially ionized gas could realize "weak reconnection dif- 
fusion". However, Lazarian et al. (2004) claimed that in the 
partially ionized gas the turbulence proceeds first as a vis- 
cosity damped magnetic turbulence with a shallow spectrum 
of magnetic perturbations and steep velocity spectrum. Then 
this regime transfers to strong turbulence involving only ions 
and electrons. From the scale of the strong turbulence res- 
urrection to the dissipation scale one expect to observe the 
regime of the accelerated diffusion making the microphysics 
on the smaller scales irrelevant. A more fundamental prob- 
lem is that in the regime of weak turbulence without a transi- 
tion to the strong one, the reconnection may get dependent on 
resistivity, i.e. become slow. In this situation one would ex- 
pect the accumulation of magnetic winding and also eventual 
transfer to bursty strong regime even in one fluid approxima- 
tion. These issues are, however, are beyond the scope of the 
present paper. For astrophysical settings that we discuss in 
§8 and §9 we do not expect the reconnection diffusion to be 
in weak regime. 
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8 Theoretical expectations and numerical simulations 
of reconnection diffusion 

8.1 Limitations of numerical simulations 

Recently we performed a few numerical studies to explore 
the consequences of reconnection diffusion. Similarly, as 
numerical studies of ambipolar diffusion do not "prove" the 
very concept of ambipolar diffusion, our studies were not in- 
tended to "prove" the idea of reconnection diffusion. Our 
goal was to demonstrate that, in agreement with the theo- 
retical expectations, the process of reconnection diffusion is 
important for a number of astrophysical set-ups relevant to 
star formation. 

We have to admit that the limitations arising from numer- 
ical simulations are not always appreciated within the as- 
trophysical community. While low resolution observations 
provide a true picture smoothed by a telescope beam, poten- 
tially, low resolution numerical simulations may provide a 
completely erroneous physical picture. 

To understand the difference between reconnection in 
astrophysical situations and in numerical simulations, one 
should recall that the dimensionless combination that con- 
trols the resistive reconnection rate is the Lundquist num- 
ber given by Eq. (l5lFl Because of the huge astrophysical 
length-scales L x involved, the astrophysical Lundquist num- 
bers are also huge, e.g. for the ISM they are about 10 16 , 
while present-day MHD simulations correspond to S < 10 4 . 
As the numerical efforts scale as L\, where L x is the size 
of the box, it is feasible neither at present nor in the foresee- 
able future to have simulations with realistically Lundquist 
numbers. Therefore our numerical studies of reconnection 
diffusion in Santos-Lima et al. (2010, 2012) deal with a 
different domain of Lundquist numbers and the theoretical 
justification why for the given problem the difference in the 
Lundquist numbers is not essential is mandatory. For the case 
of reconnection diffusion simulations, LV99 theory predicts 
that the dynamics of reconnection is independent from the 
Lundquist number and therefore the reconnection in the com- 
puter simulations in the presence of turbulence adequately 
represents the astrophysical process. 

We can rely on LV99 theory as by now its analytical pre- 
dictions have been confirmed (see Kowal et al. 2009, 2012, 
Lazarian et al. 2011, Vishniac et al. 2012; see also §4.3). 
Moreover, the correspondence of the LV99 theory and more 
recent ideas on spontaneous stochasticity developed in the 
theory of turbulence have been revealed in ELV11. These 
studies provide us with confidence that we understand what 
is happening with the magnetic field in the Santos-Lima 
(2010, 2011) simulations. 



~ The magnetic Reynolds number, which is the ratio of the mag- 
netic field decay time to the eddy turnover time, is defined using the 
injection velocity vi as a characteristic speed instead of the Alfven 
speed Va, which is taken in the Lundquist number. 
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Fig. 12. Removal of magnetic field via reconnection diffusion from 
cylindrical models of molecular clouds. In the process of simula- 
tions the density is accumulated at the center of the potential well 
(upper raw), while the magnetic field leaves the center (lower raw). 
From Santos-Lima et al. (2010). 



We must also stress that for the reconnection diffusion pro- 
cess turbulence is essential. Similarly with the Lundquist 
number issue, the perfect representation of astrophysical tur- 
bulence is not possible due to the huge difference of the 
Reynolds numbers in astrophysical fluids and in simulations. 
However, our theoretical study allows us to claim that for the 
reconnection diffusion the scales of turbulence injection or 
the scale of the transition to strong MHD turbulence (krans, 
see Eq. ( T24b ) matter the most. Therefore if the simulations 
resolve these scales one may rely on the reconnection diffu- 
sion in astrophysical settings being correctly represented. In 
terms of simulating reconnection diffusion, it is important to 
keep in mind that the LV99 model predicts that the largest 
eddies are the most important for providing outflow in the 
reconnection zone and therefore the reconnection will not be 
appreciably changed if turbulence does not have an extended 
inertial range. In addition, LV99 predicts that the effects of 
anomalous resistivity, including that arising from the finite 
numerical grid, do not change the rate of turbulent reconnec- 
tion (see more discussion in §4.2). 

8.2 Reconnection diffusion and magnetic field removal 
from clouds 

Understanding of the nature of reconnection diffusion allows 
one to simulate process with 3D MHD codes. Some results 
of such a simulation is shown in Figure [12] We observe that 
in the absence of ambipolar diffusion, the magnetic field es- 
capes the gravitational potential, allowing the matter to be- 
come concentrated in the center. The simulations were per- 
formed both from equilibrium configurations simulating sub- 
critical clouds and in collapsing clouds, simulating supercrit- 
ical clouds. In all the cases, the efficient removal of magnetic 
flux from clouds was observed. 



A. Lazarian: Reconnection and Star Formation 



23 



It worth stressing that the application of concept of recon- 
nection diffusion is not limited to accounting for the known 
observational facts. It includes predictions, e.g. the possi- 
bility of a gravitational collapse irrespectively of the degree 
of cloud ionization. As we discuss in §9.4 reconnection dif- 
fusion provides an attractive scenario for star formation in 
special cases when the expected high ionization should make 
ambipolar diffusion prohibitively slow. 

8.3 Reconnection diffusion and circumstellar accretion 
disks 

Circumstellar disks are known to play a fundamental role 
at the late stages of star formation (see Aikawa & Nomura 
2008). Observational studies revealed that embedded mag- 
netic fields in molecular cloud cores are high enough to in- 
hibit the formation of rotationally supported disks. Ambipo- 
lar diffusion is not powerful enough to induce the removal 
of magnetic fields fast enough. This motivated Shu et al. 
(2006) to talk about effects of enhanced resistivity that can 
explain the observational data (see also more recent elabora- 
tions of this idea of microscopic resistivity in Krasnopolsky 
et al. 2010 and Li et al. 2011). On the contrary, appealing 
to fast reconnection in LV99, Lazarian & Vishniac (2009) ar- 
gued that the removal of magnetic field is due to reconnection 
diffusion. 

Figure Qj] shows results of recent simulations by Santos- 
Lima et al. (2012) which, indeed, support the notion that re- 
connection diffusion is the process that is responsible for the 
removal of magnetic fields from accretion disks. The turbu- 
lence injection in the simulations was done to mimic the tur- 
bulence associated with the process of disk formation. The 
simulations also testify that without turbulence for realistic 
parameters of resistivity the formation of disks is suppressed. 
This work established the correspondence of the properties 
of the disks produced by reconnection diffusion to observa- 
tions. While additional processes may also be important for 
solving "magnetic braking catastrophe" (see Seifried et al. 
2012), the reconnection diffusion is the process that is defi- 
nitely present in turbulent fluids and therefore is an essen- 
tial part of circumstellar accretion disks models. Apart from 
being important for the dynamics of accretion disks, the in- 
tensity of turbulent magnetic field is important for the growth 
of dust particles in accretion disks (Steinacket et al. 2010) 
through the coagulation and shuttering induced by MHD tur- 
bulence (see Lazarian & Yan 2002, Yan & Lazarian 2003, 
Hoang, Lazarian & Schlickeiser 2012). Thus the process 
should be taken into account for any self-consistent model 
of circumstellar accretion disks (e.g. see Henning & Meeus 
2011), including accretion disks around massive stars (see 
Keto& Zhang 2011). 
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Fig. 13. Formation of circumstellar disks (from left to right): in hy- 
dro simulations, MHD simulations without turbulence, MHD sim- 
ulations with unrealistically high resistivity and MHD simulations 
with turbulence at the start of simulations. Reconnection diffusion 
produces realistic disks. From Santos-Lima et al. (2012). 

9 Predictions and tests for reconnection diffusion 

Reconnection diffusion is a physical process very different 
from the ambipolar diffusion. Thus it is not surprising that 
the star formation controlled by reconnection diffusion is 
very different from the picture based on the ambipolar diffu- 
sion concept. Below we outline a few predictions that follow 
from the reconnection diffusion model. 

9.1 Reconnection diffusion in interstellar diffuse gas 

A naive picture of frozen in magnetized plasma suggests that 
fluctuations of magnetic field and density should be corre- 
lated. The correlations may not be perfect, as motions along 
magnetic field lines that compress only gas are also present, 
e.g. slow modes in the media with magnetic pressure larger 
than the pressure of the ionized gas (see Cho & Lazarian 
2002, Passot & Vazquez-Semadeni 2003). 

Observations by Troland & Heiles (1986) demonstrated a 
poor correlation between the magnetic field strength and den- 
sity, which was rather unexpected at the time the work was 
done. Note that the degree of ionization of the diffuse media 
is sufficiently large to make the effects of ambipolar diffu- 
sion negligible. We discuss the approaches to this problem 
based on the works by Passot & Vazquez-Semadeni (2003) 
and Heitsch et al. (2004) in §1 1.6 and §1 1.4, respectively. 

In view of our earlier discussion in this paper, the above 
result is expected. Indeed, we discussed that the reconnec- 
tion diffusion tends to make the magnetic energy density uni- 
formly distributed in the volume. In other words, in the pres- 
ence of reconnection diffusion mixing of density fluctuations 
by turbulent eddies is present and this should mitigate any 
density-magnetic field correlations arising from simultane- 
ous compression of magnetic field and conducting gaj^l 



In comparison the models with enhanced resistivity operate 
with the resistivities that are not motivated by the known physics. 



For highly supersonic turbulence shock formation may some- 
what alter the picture and the correlation of magnetic field and den- 
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Applying our results in §5.1 one can conclude that the dif- 
fusion of magnetic field for the diffuse interstellar media is 
determined by motions at the large scale which for super- or 
transAlfvenic driving provides the coefficient of reconnec- 
tion diffusivity LVl if the scale of interest a is larger than 
the turbulence injection scale L. Similarly, for a< L the dif- 
fusion coefficient is ~ aV^iaj 'L) 1 / 3 if a < L. For the sub- 
Alfvenic driving, i.e. for the Alfven Mach number of tur- 
bulence Ma < 1 additional factor of M\ is expected in the 
diffusion coefficient expressions reflecting the decreased ef- 
ficiency of mixing by subAlfvenic turbulence. The results of 
numerical simulations are consistent with these expectations 
(see Santos-Lima et al. 2010). In fact, reconnection diffusion 
was demonstrated to decorrelate magnetic field and density, 
if initially this correlation was present. 

Reconnection diffusion can explain other observations as 
well. For instance, Crutcher (2012) analyzed an extensive 
set of published as well as unpublished Zeeman measure- 
ments and showed that clouds with column densities Njj less 
than 10 21 cm -2 are subcritical, while at higher densities they 
get supercritical. He noted, that that for cold HI clouds with 
Nh < 10 21 cm~ 2 the magnetic field is approximately 6/iG, 
which corresponds to the value of magnetic field strength in 
a much more rarefied warm neutral media. From this he con- 
cluded that the diffuse clouds had to be formed either via the 
compression along magnetic lines or that the formation of 
clouds proceeds selectively only at the regions of low mag- 
netic field strength. In contrast to this, on the basis of our 
study we claim that reconnection diffusion presents a more 
appealing alternative. Indeed, as we discussed above, the 
efficient diffusion of magnetic field from the clouds should 
make magnetic field strength the same in cold dense clouds 
and surrounding warm rarefied medium. 

In addition, Crutcher (2012) notes that the kinetic and 
magnetic energies of such low density clouds are approxi- 
mately in equipartition and larger that the thermal energy. 
From the point of view of our analysis, this is a signature of 
the transAlfvenic supersonic turbulence indicating that the 
turbulence should be efficient in driving reconnection diffu- 
sion. From the plot presented in Crutcher (2012), namely, 
figure 7 in his paper, it is evident that the majority of the 
clouds in his compilation still preserves the same magnetic 
field strength of the order of 6/iG even as the column den- 
sity gets as high as 10 23 cm" 2 . We believe that this is the 
consequence of the reconnection diffusion being efficient at 
those densities. We see a clear tendency of the increase of 
the mean magnetic field in the sample for densities larger 
than 10 23 cm -2 . A possible explanation for this is that for 
those clouds self-gravity gets important and therefore the re- 
connection diffusion fails to remove magnetic field from the 



sity is being observed in supersonic simulations (Burkhart et al. 
2009). However, the bulk of the turbulence in warm diffuse media 
for which the observations are applicable is subsonic (see Burkhart 
et al. 2012). 



contracting clouds fast enough. Note, that the high density 
clouds tend to be smaller and therefore, as we discussed in 
§5.2 the turbulent scales involved in the reconnection diffu- 
sion get also smaller. The reconnection diffusion for such 
clouds slows down. 

We note that the observational results in Crutcher (2012) 
on the independence of magnetic field strength on the cloud 
density (and the reconnection diffusion concept that accounts 
for these results) can be related to the empirical Larson rela- 
tions (Larson 1981) obtained for interstellar turbulence. Lar- 
son (1981) found that the velocity dispersion is proportional 
to the square root of the cloud size, i.e. cry ~ i? 1 / 2 and that 
the 3D density of cloud is inversely proportional to cloud 
size, i.e. p^ R~ l , For instance, one can assume a rough 
equality between the kinetic energy and magnetic energy 
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which is natural for transAlfvenic turbulence as well as the 
virialization of a cloud 
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Combining Eqs. d32b and ( 1331 with a simplest estimate of 
cloud mass M ~ pR 3 one gets 



and 
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which reproduce the Larson (1981) relations if the reconnec- 
tion diffusion keeps magnetic field uncorrelated with den- 
sity. For cloud cores where reconnection diffusion is not fast 
enough to remove magnetic field on the time of the dynamic 
collapse (see Tafalla et al. 1998, Reiter et al. 201 1), the Lar- 
son relations fail, in agreement with observations and simula- 
tions (see Nakamura & Li 201 1). Quantitatively, the criterion 
for the reconnection diffusion to be able to remove the mag- 
netic field from the collapsing cloud is Vi n f a ii < k/I c IoucI, 
where k is one of the diffusion coefficients the choice of 
which depends on the regime of the turbulence (see §5.1). 

9.2 Core and envelope magnetization 

Recently the idea of star formation being mediated by am- 
bipolar diffusion was challenged in Crutcher et al. (2010, 
2011). The authors measured matter magnetization in the 
core of a dark cloud and in the envelope surrounding the 
cloud (see our schematic representation in Figure [T4l). Con- 
trary to the predictions of the theory based on ambipolar 
diffusion, the observations showed that the core is more 
magnetized that the envelope. Naturally, this work, which 
challenged the key predictions of the ambipolar diffusion 
paradigm induced intensive controversy (Mouschovias & 
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Tassis 2009, 2010). Without getting into details of the well- 
publicized arguments of both sides, let us pose a question of 
whether the results observed by Crutcher et al. (2010, 201 1) 
are consistent in the presence of reconnection diffusion. 

As we mentioned earlier, in the presence of gravity recon- 
nection diffusion allows density to concentrate towards the 
gravitational center, while magnetic field leave the gravita- 
tional potential. The reconnection diffusion rate and thus the 
rate of magnetic field removal is expected to be proportional 
to the level of turbulence (see Eqs. (fT8l l. (l25l l. (l26li). Simula- 
tions that illustrate this effect are presented in Figure[l4] It is 
known that the velocities in cloud envelopes are larger than 
the velocities in the cloud cores (see Taffalla et al. 1998). In 
addition, the scale of turbulence involved in reconnection dif- 
fusion is also larger. Thus we expect the diffusion coefficient 
for the envelope (see Eq. l25l l to be larger for the envelope 
compared to the core and faster removal of magnetic field 
from the envelope than from the core. This agrees with the 
results in Crutcher et al. (2010). 

9.3 Predictions for the big picture of star formation 

Both reconnection diffusion and ambipolar diffusion remove 
magnetic field and can initiate gravitational collapse of the 
gas. However, there are substantial differences between the 
two processes. 

The characteristic diffusion time for magnetic fields is 
given by K m /L 2 , where n m is given by Eqs. ( TT8b , d25l >. d26b 
depending on the regime of turbulence. The diffusion coef- 
ficient K m varies spatially, as the turbulence changes from 
place to place and also with density. For individual clouds 
and cores if we neglect, for the sake of simplicity, the ef- 
fects of stratification the reconnection diffusion becomes the 
function of turbulent driving. For subAlfvenic driving one 
gets a suppression factor of M\, which reflects the ineffi- 
ciency of weak turbulence to induce magnetic field diffusion. 
For superAlfvenic and transAlfvenic driving, the magnetic 
field is to be transported with the turbulent diffusivity rate. 
The length scale of turbulence depends on the sources of the 
turbulence. If the turbulence is driven externally, the corre- 
sponding length-scale is expected to be of the order of the 
cloud size, if the stratification of the cloud is negligible. 

In diffuse interstellar gas and Giant Molecular Cloud 
(GMC) the reconnection diffusion is expected to dominate. 
One may provide evidence in support of the reconnection 
diffusion concept. For instance, unlike ambipolar diffusion, 
the reconnection does not depend on atomic level or dust 
physics. It was pointed to us by Bruce Elmegreen (private 
communication) that it is observed that the star formation is 
about the same in galaxies with low metallicities as in galax- 
ies with high metallicities, which is hard to understand if 
magnetic field loss is governed by ambipolar diffusion. In- 
deed, the latter is supposed to be much faster in low-Z galax- 
ies with high metallicities and, within the ambipolar diffu- 
sion paradigm it is expected, contrary to observations, that 



the Initial Mass Function (IMF) and star formation rate to 
shifted in low-Z galaxies. This is not a problem for the re- 
connection diffusion which is controlled only by the media 
turbulence. 

Within the standard star formation paradigm the low ef- 
ficiency of star formation is bottlenecked by the ambipolar 
diffusion rate. This paradigm has been altered recently and 
the low efficiency of star formation was attributed to the low 
fraction of the cloud mass that has sufficiently high density to 
form stars (Elmegreen 2007, Vazquez-Semadeni et al. 2009). 
In the latter picture, star formation blows the cloud cores 
apart, while the efficiency of star formation is high in the 
cloud core with up to third of the mass going into stars. At 
the same time, in terms of the mass of GMCs the efficiency 
is low e.g. less than five percent. In other words, the cloud 
envelope is passive in terms of star formation, the matter is 
being dispersed or/and pushed aside without forming stars. 
In the GMC envelopes, the reconnection diffusion is going 
to be efficient making sure that density and magnetic fields 
are well mixed. On the available time scales it will tend to 
equalize the magnetic field strength within the clouds and 
the ambient interstellar medium, inducing the loss of part of 
the flux captured in the cloud at the stage of its formation. 
An important difference with the ambipolar diffusion is that 
the reconnection diffusion is not only associated with the re- 
moval of the magnetic field, but also with its turbulent mix- 
ing which tends to make the distribution of magnetic field 
uniform. Turbulent mixing helps keeping the diffuse media 
in magnetized subcritical state. In this situation, the exter- 
nal pressure is important for initiating collapse and which 
well corresponds to the observations of numerous small dark 
clouds not forming stars in the inter-arm regions of galax- 
ies (Elmegreen 201 1). The reconnection diffusion process is 
expected to relax the sharp local changes of magnetic field 
direction, thus explaining the alignment of magnetic field of 
molecular clouds cores with the magnetic field of the spiral 
arms (Li & Henning 201 1). 

Ambipolar diffusion depends on the ionization of mate- 
rial and this results in the introduction of characteristic den- 
sity for the clouds and cores to become leaky to magnetic 
field. In contrast, the reconnection diffusion concept im- 
plies that in realistic turbulent media there is no character- 
istic density for the collapse to be initiated. Therefore any 
cloud with the appropriate virial parameter (see McKee & 
Zweibel 1992) can form stars. The difference between dif- 
ferent clouds stems from the density controlling the timescale 
of the collapse, level of steering turbulence or temperature of 
media. Directly, the requirement of clouds to be molecular is 
not present for the reconnection diffusion to induce star for- 
mation. However, molecular clouds have lower temperature. 

We would like to emphasize that the role of turbulence in 
the presence of reconnection diffusion is two fold. First of 
all, it allows the removal of magnetic flux, stimulating star 
formation within a contracting cloud. However, solenoidal 
steering also changes the virial balance making clouds less 
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Fig. 14. it Left panel. Schematic of the cloud core and envelope in Crutcher et al. (2010). Both scale of the turbulence and velocity (and thus 
reconnection diffusion rate) are larger in the envelope than in the core. Right panel. Results of simulations (modified from Santos-Lima et 
al. 2010) The removal of magnetic field via reconnection diffusion increases with the increase of the turbulence level. 



prone to collapse. The amplitude of the steering motions in- 
creases as vi ~ Z 1 / 3 making larger, e.g. diffuse atomic clouds, 
not eligible for a gravitational collapse. The compressible 
motions associated with turbulence at the same time increase 
the density and may stimulate the collapse. The ratio be- 
tween compressible and incompressible motions depends on 
turbulence driving. However, simulations (e.g. Cho & Lazar- 
ian 2002) testify that the most of energy tend to reside in 
solenoidal motions. Therefore, one would expect that the 
increase of the level of turbulence decreases star formation, 
although it cannot shut it down completely. A further discus- 
sion of this issue is beyond the scope of this paper. It is clear 
that reconnection diffusion reveals a new important role that 
is played by interstellar turbulence. 

9.4 Reconnection diffusion and extreme cases of star 
formation 

Reconnection diffusion concept provides new ways of ap- 
proaching challenging physics of star formation in extreme 
environments. For instance, galaxies emitting more than 
10 12 solar luminosities in the far-infrared are called ultra- 
luminous infrared galaxies or ULIRGs. The physical con- 
ditions in such galaxies are extreme with very high density 
of cosmic rays (see Papadopoulos et al. 2011). Ambipo- 
lar diffusion is expected to be suppressed due to cosmic ray 
ionization. At the same time these environments provide the 
highest star formation rate which is suggestive of a process 
which removes magnetic fields irrespectively of the level of 
ionization. Reconnection diffusion is such a process. 

The formation of early stars (see Chiappini et al. 2011) 
is a great problem for which the effects of magnetic fields 
are debated. Reconnection diffusion mitigates the effects of 



magnetic fields and therefore decreases the uncertainties as- 
sociated with the presence of magnetic fields at the sites of 
early Universe star formation. 

10 Additional consequences of reconnection diffusion 

10.1 Reconnection diffusion and cosmic ray accelera- 
tion in dark clouds 

The lower energy e.g. ~ 100 MeV cosmic rays dominate the 
ionization of cool neutral gas, especially dark UV shielded 
molecular regions (Goldsmith & Langer 1978). Some obser- 
vations indicate substantial changes of the cosmic rate ion- 
ization rate between diffuse and dense gas (see McCall et al. 
2003) which may be the consequence of local cosmic ray ac- 
celeration. As we discussed, reconnection diffusion invokes 
LV99 model of reconnection. However, the same type of 
reconnection is expected to accelerate cosmic rays (de Gou- 
veia dal Pino & Lazarian 2005, Lazarian 20050 The ac- 
celeration happens as particles bounce back and forth within 
shrinking magnetic loops (see Figure \15[. Recently, the ac- 
celeration of cosmic rays in reconnection has been invoked 
to explain results on the anomalous cosmic rays obtained by 
Voyager spacecrats (Lazarian & Opher 2009, Drake et al. 
2010), the local anisotropy of cosmic rays (Lazarian & De- 
siati 2010) and the acceleration of cosmic rays in clusters of 
galaxies (Lazarian & Brunetti 2011). Naturally, the process 
of acceleration is much wider spread and not limited by the 
explored examples. 



The predicted spectrum without taking the backreaction of the 
accelerated particles is N(E)dE ~ E^^^dE. Considerations in 
Drake et al. (2006) suggest that the spectrum of the particles can 



A. Lazarian: Reconnection and Star Formation 



27 



v R t 




v4 

Fig. 15. First order Fermi acceleration as cosmic rays bounce within 
a 3D loop of reconnected flux that shrinks due to magnetic recon- 
nection. From Lazarian 2005. 



The maximal energy of the particles can be estimated 
with the standard formulae (see Longair 2011) E max ~ 
10 16 eVBR, where magnetic field B is normalized by 1 /xG 
and the radus is normalized by 1 pc. This maximal energy 
may not be reached due to the CR losses if the acceleration 
is insufficiently efficient. Numerical simulations by Kowal et 
al. (2012) show a relative inefficiency of first order Fermi ac- 
celeration by reconnection in pure turbulence without a large 
scale magnetic reversal. This is likely to be due to the gen- 
tle variations of magnetic field in the GS95 model of tur- 
bulence. As we discussed earlier (see §5.1) reconnection in 
strong MHD turbulence does necessarily take place, but the 
reconnecting bundles of magnetic field cross each other at a 
small angle </> and the resulting rate of energy gain Va cos<j>/c 
is low. Gravitational forces acting on a collapsing core suffi- 
ciently perturb magnetic field lines increasing the efficiency 
of the acceleration within reconnection regions. 

Stretching of the magnetic field lines happens both inside 
molecular clouds and outside them as curved magnetic field 
lines diffuse away. The Alfven velocity outside clouds is 
larger and this may provide a more efficient acceleration in 
agreement with the ionization results for the ( Persei cloud 
(see Le Petit 2004). Naturally, more 3D modeling^ of the 
acceleration in reconnection regions is necessary to seek the 
quantitative agreement with observations. 



get shallower if the backreaction is taken into account. 

31 Modeling in Kowal et al. (201 lb) showed that the acceleration 
in 2D and 3D proceed at a different rate, which questions the ap- 
plicability of results obtained in 2D simulations (see Drake et al. 
2010) to the actual astrophysical systems. 



10.2 Heating of clouds 

Ambipolar drag induces additional heating of star forming 
clouds, which can potentially be detected as the excess of the 
heat after turbulent, cosmic ray etc. ways of heating are ac- 
counted for. As any fast reconnection model, LV99 predicts 
that only a small fraction of magnetic energy is going to be 
dissipated directly through Ohmic heating. Most of the en- 
ergy is going to be consumed by the turbulent outflow. Thus 
in the process of reconnection diffusion we expect that the 
heating is going to be linked with turbulence in clouds. Ad- 
ditional turbulent energy input might be detected through the 
variations of the power spectrum. The latter can be measured 
by the Velocity Channel Analysis (VCA) or Velocity Cor- 
relation Spectrum (VCS) techniques (Lazarian & Pogosyan 
2000, 2004, 2006, Padoan et al. 2004, 2009, Chepurnov & 
Lazarian 2009, Chepurnov et al. 2010, see Lazarian 2010 for 
a review). 

10.3 Star formation simulations: criterion for repre- 
senting magnetic diffusion 

We believe that flux freezing being the dominant idea in as- 
trophysical community strongly influenced the interpretation 
of the results of numerical simulations. Indeed, within the 
flux freezing paradigm the only way to explain how conduc- 
tive matter become concentrated without dragging magnetic 
field lines along is to postulate its motion along magnetic 
field lines. Any other explanation would open a a pandora 
box of the necessity to justify how diffusive numerical sim- 
ulations can represent highly conductive astrophysical envi- 
ronments. 

The concept of reconnection diffusion provides a new out- 
look at the existing simulations. Plasmas and magnetic field 
are subject to efficient reconnection diffusion both in astro- 
physical conditions and in numerical simulations. As we dis- 
cussed in §4.2 and §6.4 the rate of reconnection diffusion is 
independent of the detailed microphysics of local reconnec- 
tion events. In this respect, numerical small scale reconnec- 
tion does not compromise the results of numerical simula- 
tions of turbulent environments. However, at scales where 
numerical diffusion suppresses turbulence the physics repro- 
duced by codes is different in simulations and in real astro- 
physical conditions. 

The field wandering weakly which is essential for the 
LV99 reconnection depends on the amplitude and scale of 
turbulent motions and does not depend on the small scale 
microphysics. These is excellent news for MHD simulations 
of astrophysical turbulent environments, in particularly, the 
star formation. As we mentioned in §8.1 astrophysical simu- 
lations are not in the right ballpark of the magnetic Reynolds 
or Lundquist numbers and it is important that in terms of 
magnetic reconnection this huge difference do not alter the 
physics of the process. 
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In fact, reconnection diffusion has always been part of the 
numerical simulations involving turbulence even if its role 
has not been identified. For instance, Julian Krolik pointed 
to us his work on magnetic field transport around a black hole 
(Beckwith, Hawley & Krolik 2009). We also believe that re- 
connection diffusion is the reason for the expulsion of the 
magnetic flux in the simulations of Romanova et al. (201 1). 
We might speculate the simulations in Li et al. (2012) are 
also affected by reconnection diffusion. We may argue that 
some effects explained as the consequence of the accumu- 
lation of matter along magnetic field lines (see §11.6) are 
also due to reconnection diffusion. However, without de- 
tailed studies of the aforementioned simulations it is not 
clear whether the real reconnection diffusion physics or bo- 
gus effects of numerical diffusivity are at play. For instance, 
Crutcher et al. (201 1) refers to the simulations in Luntilla et 
al. (2009) which produce, in agreement with observations, 
higher magnetization of the cloud cores. If these cores are of 
the size of several grid units across, numerical effects rather 
than reconnection diffusion may be dominant and turbulence 
is suppressed at these scales. 

All in all, one can formulate the criterion for the numeri- 
cal simulations to represent the actual diffusion of magnetic 
fields: if on the scales of study simulations exhibit turbulence, 
it is reconnection diffusion that dominates the numerical one. 
If this criterion is satisfied, in terms of the magnetic field dif- 
fusion, the simulations are trustworthy. Thus we claim that 
if the cores in the simulations are turbulent, the removal of 
magnetic flux from them via reconnection diffusion is similar 
to astrophysical highly conductive environment. This opens 
a possibility of developing adaptive mesh codes designed to 
reproduce astrophysical magnetic diffusion correctly in spite 
of the difference in Lundquist numbers. 

11 Discussion 

11.1 Physics of reconnection diffusion 

The idea that fast reconnection predicted by LV99 model can 
change the picture of star formation by inducing magnetic 
field diffusion has been discussed in a number of papers start- 
ing with Lazarian (2005). However, the physical picture of 
the process has not been properly discussed. This was a prob- 
lem, as the competing picture of magnetic drift via ambipolar 
diffusion is very clear in terms of the physics involved (see 
Mouschovias 1991). 

The above shortcoming is dealt with in the present paper 
where two physical pictures of the process, macrophysical 
based on LV99, and microphysical, based on ELV11, are 
presented. The two approachers are equivalent in the very 
core of the underlying physics, but in terms of understanding 
of the process they exhibit complementary features. Indeed, 
LV99 primary deals with large scale fluxes that intersect and 
reconnect. This is the standard way of dealing with magnetic 



reconnection. The large scale variations of magnetic field 
must be dealt in star formation and the answer provided by 
LV99 is that such variations do not stop the process of recon- 
nection diffusion. Naturally, as we also showed in this paper, 
the LV99 model can be applied to different scales of mag- 
netic field hierarchy, including the hierarchy of the turbulent 
eddies within the magnetized flow. At this point it enters the 
domain of self consistently describing the cascade of recon- 
nection events in the turbulent flow that is important for the 
continues changes of magnetic field line identities that is es- 
sential for reconnection diffusion. This is also the domain 
which is described by the concept of "spontaneous stochas- 
ticity" (see ELV1 1). The latter, instead of describing the con- 
stantly and stochastically changing magnetic field line con- 
nections, uses the description which explicitly appeals to the 
magnetic field stochasticity at turbulent scales and the viola- 
tion of flux freezing. 

The gain in physical understanding of the reconnection 
diffusion enabled us to provide a better description of the 
consequences of the process of star formation. For instance, 
it allowed us to better identify the role of weak and strong tur- 
bulence for transporting magnetic fields. The estimates that 
we provided can be tested both observationally and numeri- 
cally. In view of the latter point, we would like to stress that 
by itself the existing numerical studies of reconnection diffu- 
sion cannot be used to prove the concept. The simulations so 
far have been done in MHD regime and instead of the gener- 
alized Ohm's law (see Eq. (l28l i have used ordinary resistiv- 
ity. In this situation it is the LV99 reconnection theory and 
the physically equivalent to it theory of spontaneous stochas- 
ticity (ELV11) that justify that why one can disregards the 
small scale physics dealing with the reconnection diffusion 
process. 

Our results in the paper should not be understood as the 
claim that ambipolar diffusion is unimportant for star for- 
mation. Future research, both observational and theoreti- 
cal, should determine the relative importance of reconnection 
and ambipolar diffusion processes. As observations show 
that turbulence is ubiquitous in diffuse ISM and molecular 
clouds, we expect the reconnection diffusion to be important 
for magnetic field transport in these environments. However, 
in very quiescent cloud cores ambipolar diffusion may be the 
dominant process. 

11.2 Intuitive understanding of reconnection diffusion 

The idea of magnetic flux freezing is so deeply rooted in 
astrophysics, that any attempts to challenge it sound rather 
heretical. We would like to stress that fast reconnection in 
turbulent magnetized media is very natural. Without it the 
turbulent fluid would create felt-like structures, as was advo- 
cated by Don Cox (private communication). Indeed, inter- 
secting and not being able to pass through each other mag- 
netic field lines are bound to stop magnetized fluid behave 
like a fluid. If, however, this is not the case and magne- 
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tized fluids preserve fluid-type behavior in the presence of 
turbulence, one has to accept the efficient diffusion of mat- 
ter and magnetic fields. This happens through reconnection 
diffusion as we explained in this paper. It worth mentioning 
that the ideas of magnetic field meandering that have been 
invoked for decades to understand the observed diffusion of 
cosmic rays perpendicular to magnetic field (see Parker 1965, 
Jokipii 1973) are naturally related to the reconnection diffu- 
sion concept. 

As we discussed in the paper, the reconnection diffusion is 
closely related to the flux freezing violation in turbulent envi- 
ronments. The failure of flux freezing in turbulent fluids has 
very big astrophysical consequences for star formation and 
beyond it. The widely accepted point of view is that in astro- 
physical situations the flux freezing is "nearly" fulfilled and 
the violations are due to the existing finite non-ideal effects. 
The problem is that this is not true in realistically turbulent 
astrophysical fluids. However, one should keep in mind that 
turbulence, unlike resistivity, does not destroy magnetic field 
lines, but it makes the magnetic field stochastic. The total 
magnetic flux does not change, but the charged particles get 
the possibility of exploring the entire volume, which also 
means that magnetic field gets diffusive. In turn, the latter 
entails magnetic field filling the entire volume and mitigat- 
ing its effects of counteracting gravitational compressing the 
matter. 

11.3 Turbulent and magnetic support of clouds 

The traditional approach to star formation frequently appeals 
to turbulent and magnetic support of the clouds (see Mestel 
1965). As magnetic field is usually assumed to be in rough 
equipartition with kinetic motions, one may wonder whether 
neglecting magnetic support just amounts to the factor of or- 
der unity in the picture above. 

There is, however, a serious difference between the effects 
of magnetic field and turbulence. For the virial support by 
turbulence, only motions less than the cloud scale are impor- 
tant. The larger scale motions do not enter the virial equation. 
The outside cascade, if anything, can compress the cloud due 
to compressible turbulent fluctuations and cloud-turbulent in- 
teractions. The role of magnetic field is different. If the 
large scale magnetic field is dragged into the cloud, it can 
only be amplified by the cloud compression. Therefore the 
large scale equipartition between turbulence and magnetic 
field does not preclude the magnetic field to be dominant on 
the scale of self-gravitating cores. 

The process of reconnection diffusion allows the magnetic 
field to equalize inside and outside the cloud, decreasing the 
effect of magnetic support. This corresponds to the modern 
understanding of star formation as a very dynamic process 
with no necessity of support of small scale infall of matter 
(see Elmegreen 201 1 and references therein). 

Incidentally, reconnection diffusion shows that the text- 
book picture of the magnetized cloud with magnetic field 



support in the direction perpendicular to magnetic field and 
turbulence providing the vertical extend of the cloud is not 
sustainable. Reconnection diffusion is expected to remove 
the excess of magnetic field from the cloud on the dynamical 
time scales. The conventional picture above holds only in the 
absence of strong MHD turbulence. 

11.4 Reconnection diffusion and turbulent ambipolar 
diffusion 

An interesting study focused on the ambipolar diffusion 
physics in a turbulent flow was performed Heitsch et al. 
(2004, henceforth HX04). They performed 2.5D simulations 
of turbulence with two-fluid code and examined the decorre- 
lation of neutrals and magnetic field that was taking place as 
they were driving the turbulence. The study reported an en- 
hancement of ambipolar diffusion rate compared to the am- 
bipolar diffusion acting in a laminar fluid. HX04 correctly 
associated the enhancement with turbulence creating density 
gradients that are being dissolved by ambipolar diffusion (see 
also Zweibel 2002). Due to magnetic field being perpen- 
dicular to the flow in 2.5D, the set up precluded magnetic 
field reconnection and magnetic fields preserved their iden- 
tity throughout the simulations (cf. §6). Thus HX04 studied 
the effect different from reconnection diffusion that we deal 
in the paper. They termed the process "turbulent ambipolar 
diffusion". 

The set up in HX04 presents a special case of a magne- 
tized flow when magnetic fields act only as an additional 
pressure within the fluid that exhibits hydrodynamic behav- 
ior irrespectively of the strength of magnetic field. Therefore, 
similar to the hydrodynamic case, turbulence in this situation 
cannot be weak, even if Va ^> Vl (cf. Table 1) and 2.5D ed- 
dies stir the fluid irrespectively of the strength of the driving. 
Ions stayed entrained on the field lines (cf. Figure |9]l while 
both neutral density and magnetic field spread diffusively at 
approximately the eddy rate. The decorrelation arose due to 
small scale ion-neutral drifts. 

Both viscosity and the decorrelation of magnetic field and 
neutral density in HX04 are due to the slippage of neutrals 
and ions and therefore in the presence of hydrodynamic-type 
eddies the diffusivity given by Eq. (TTSI l agrees well with 
their finding. In analogy with hydrodynamics, we expect that 
the diffusivity would not change even as ambipolar diffusion 
rate changes. Indeed, if the rate of ambipolar diffusion gets 
smaller the turbulent cascade proceeds to smaller scales al- 
lowing mixing at those scaleo- The reduced ambipolar dif- 

32 A similar process takes place in the case of molecular diffu- 
sivity in turbulent hydrodynamic flows. The result for the latter 
flows is well known: in turbulent regime molecular diffusivity is 
irrelevant for the turbulent transport. Indeed, in the case of high 
microscopic diffusivity, the turbulence provides mixing down to a 
scale h at which the microscopic diffusivity both, suppresses the 
cascade and ensures efficient diffusivity of the contaminant. In the 
case of low microscopic diffusivity, turbulent mixing happens down 
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fusion would still be adequate to decorrelate the magnetic 
field and the reduced scales. We infer that the limiting case 
of this flow is the 2.5D flow with no ambipolar diffusion but 
with the diffusivity still given by Eq. ( fTol l. This still would 
not be the reconnection diffusion case, as no reconnection is 
allowed. 

The difference of what we discussed in the present pa- 
per and the idea in HX04 arises from the difference in the 
flows that we considered. The set up that we deal with is 
a generic 3D turbulent flow where magnetic reconnection 
is inevitable. Therefore, as we discussed in the paper, ions 
themselves are diffusive (see §6) and the ionic density decor- 
relates with magnetic field as well. Similarly to what we 
discussed earlier, in the absence of ambipolar diffusion, the 
turbulence propagates to smaller scales making small-scale 
interactions possible. On the other hand, ambipolar diffusion 
affects the turbulence, increasing its damping. As a result, 
analogously to 2.5D case above, the ambipolar diffusion acts 
in two ways, in one to increase the small-scale diffusivity 
of the magnetic field, in another is to decrease the turbulent 
small-scale diffusivity and these effects essentially compen- 
sate each otheiF 3 ! Therefore, we believe that, at least for the 
case of strong MHD turbulence, ambipolar diffusion does not 
play any role for the turbulent transport in magnetized fluid. 
In a generic situation of 3D magnetized turbulence, recon- 
nection is essential and reconnection diffusion takes place. 
The diffusion coefficient ~ LVt corresponds to our predic- 
tion of reconnection diffusion induced by transAlfvenic and 
super Alfvenic turbulence. We note that, in the presence of 
turbulence, the independence of the gravitational collapse 
from ambipolar diffusion rate was reported in numerical sim- 
ulations by Balsara, Crutcher & Pouquet (2001). 

What does happen in the subAlfvenic case of weak turbu- 
lence? As we mentioned before, we are not aware of the stud- 
ies of diffusion in two fluids in this regime. At the same time, 
this may be potentially the most interesting case as far as the 
interplay of turbulence and ambipolar diffusion is concerned. 
As we discussed in §3.2 and §5.1 MHD turbulence at large 
scales corresponds to the "weak" regime and can be viewed 
as a collection of non-linear weakly interacting waves. How- 
ever, at a smaller scale, namely at the scale hrans given by 
Eq. d24l > the turbulence gets into the regime of strong inter- 
actions, when the intensive mixing happens in the direction 

to a scale 1% <C£i, which ensures that even low microscopic diffu- 
sivity is sufficient to provide efficient diffusion. In both cases the 
total effective diffusivity of the contaminant is given by the product 
of the turbulent injection scale and the turbulent velocity. 

33 A possible point of confusion is related to the difference of the 
physical scales involved. If one associates the scale of the reconnec- 
tion with the thickness of the Sweet-Parker layer, then, indeed, the 
ambipolar diffusion scale is much larger and therefore the recon- 
nection scale gets irrelevant. However, within the LV99 model of 
reconnection, the scale of reconnection is associated with the scale 
of magnetic field wandering. The corresponding scale depends on 
the turbulent velocity and is not small. 



perpendicular to the local direction of magnetic field. The 
diffusivity associated with turbulence is given by Eq. (1251 . 
If the ambipolar diffusivity is less than this value, it will not 
play any role and the diffusivity will be purely "turbulent". 
If, however, damping happens at scale larger than t trans a hy- 
pothetical new regime of "weak reconnection diffusion" may 
be present (see §7). The study of "turbulent ambipolar dif- 
fusion" with weak MHD turbulence has not been performed 
as far as we know. The effects of the enhancement of the 
total diffusivity are thus unclear. We might expect little, if 
any, parameter space for the "turbulent ambipolar diffusion" 
when turbulence and ambipolar diffusion synergetically en- 
hance diffusivity, acting in unison. Nevertheless, this point 
should be tested by three-dimensional two-fluid simulations 
exhibiting both ambipolar diffusion and turbulence. The ef- 
fect to search in order to test the effect of "turbulent ambipo- 
lar diffusion" is the enhancement of the diffusivity compared 
to the rate of diffusivity arising from the weak subAlfvenic 
turbulence (given by Eq. d23l). 

11.5 Reconnection diffusion and hyper- resistivity 

To explain fast removal of magnetic field from accretion 
disks Shu et al. (2006) appealed to the hyperrestivity concept 
(Strauss 1986, Bhattacharjee & Hameiri 1986, Hameiri & 
Bhattacharjee, Diamond & Malkov 2003). The studies intro- 
ducing hyperresistivity attempt to derive the effective resis- 
tivity of the turbulent media in the context of the mean-field 
resistive MHD. Using magnetic helicity conservation the au- 
thors derived the electric field. Then, integrating by part, 
one obtained a term which could be identified with effec- 
tive resistivity proportional to the magnetic helicity current. 
There are several problems with this derivation. In particular, 
the most serious is the assumption that the helicity of mag- 
netic field and the small scale turbulent fields are separately 
conserved, which erroneously disregard the magnetic helic- 
ity fluxes through open boundaries, which is essential for fast 
stationary reconnection (see more discussion in Kowal et al. 
2009 and ELV11). In more general terms, hyper-resistivity 
idea is an incarnation the mean-field approach to explaining 
fast reconnection. As explained in ELV11, the problem of 
such approaches is that the lines of the actual astrophysical 
magnetic field should reconnect, not the lines of the mean 
field. Therefore the correct approach to fast reconnection 
should be independent of the spatial and time averaging. 

All in all, we believe that the concept of hyperresistivity 
is poorly justified and should not be applied to astrophysical 
environments. 

1 1 .6 Reconnection diffusion and collecting matter along 
magnetic field 

As we mentioned in §9.1 an alternative way for changing 
flux to mass ratio is to allow conducting matter to be accumu- 
lated along magnetic field lines. This process definitely takes 
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place, but the prescription of one dimensional motion of mat- 
ter is very restrictive. More importantly, as we discussed in 
§6 the idea of fixed magnetic field lines is not applicable to 
turbulent magnetized fluid. Thus the effects of reconnec- 
tion diffusion should inevitably interfere even if plasmas is 
launched along magnetic field lines. 



For instance, our work on reconnection diffusion in diffuse 
interstellar medium should also be distinguished from the re- 
search on the de-correlation of magnetic field and density 
within compressible turbulent fluctuations. Cho & Lazar- 
ian (2002, 2003) performed three-dimensional MHD simu- 
lations and reported the existence of separate turbulent cas- 
cades of Alfven and fast modes in strongly driven turbulence 
as well as a cascade of slow modes driven by Alfvenic cas- 
cade. Slow modes in magnetically dominated plasma are 
associated with density perturbations with marginal pertur- 
bation of magnetic fields, while the same is true for fast 
modes in weakly magnetized or high beta plasmas. Natu- 
rally, these two modes de-correlate magnetic fields and den- 
sity on the crossing time of the wave. This was the effect 
studied in more detail in one-dimensional setting both both 
analytically and numerically by Passot & Vazquez-Semadeni 
(2003), who stressed that the enhancements of magnetic field 
strength and density may correlate and anti-correlate in tur- 
bulent interstellar gas within the fluctuations and this can in- 
troduce the dispersion of the mass-to-flux ratios within the 
turbulent volume. Each of the fluctuations provide a tran- 
sient change of the pointwise magnetization. In the absence 
of other effects, e.g. related to the thermal instability, the 
de-correlation is reversible. In comparison, the "reconnec- 
tion diffusion" deals with the permanent de-correlation of 
magnetic field and density making magnetic field-density de- 
correlation irreversible. In many instances both processes act 
together providing the observed (see §9.1) decorrelated mag- 
netic field and density state. 



While for the diffuse media both collecting the matter 
through compressible modes and mixing through reconnec- 
tion diffusion can act together, the idea of star formation 
based on collecting matter along magnetic field lines is more 
problematic. The scales for such one-dimensional collection 
are enormous (~ 1 kpc see Vazquez-Semadeni et al. 2011) 
and it is not feasible that in turbulent environments reconnec- 
tion diffusion would not interfere with the postulated one- 
dimensional motion. We believe that reconnection diffusion 
is an intrinsic part of the simulations of the turbulent inter- 
stellar medium and it should be accounted for in interpreting 
the results of numerical simulations provided that the crite- 
rion in §10.3 is satisfied. If the criterion is not satisfied, then 
the magnetic diffusivity is dominated by numerical effects 
and one should be cautious interpreting such simulations. 



11.7 Reconnection diffusion and modern understand- 
ing of MHD turbulence and reconnection 

The concept of reconnection diffusion is deeply rooted in the 
modern understanding of MHD turbulence and its intrinsic 
connection with reconnection. While mixing motions per- 
pendicular to the local direction of magnetic field is a part of 
the GS95 picture of Alfvenic turbulence, the star formation 
textbooks frequently depict magnetic turbulence as a collec- 
tion of waves with wavevectors mostly parallel to magnetic 
field. The latter induce only marginal mixing of matter and 
therefore other mechanisms of changing flux to mass ratio 
are required. 

Mixing motions inevitably induce the question of mag- 
netic reconnection. Fortunately, as we discussed at the end 
of section §4.3 the LV99 model provides the necessary rates 
thus providing a physical justification for both the GS95 
model and the reconnection diffusion concept. 

In the paper we considered MHD turbulence where the 
flows of energy in the opposite directions are balanced. 
When this is not true, i.e. when the turbulence has non-zero 
cross-helicity, or imbalanced (see §4.3). Studies of recon- 
nection diffusion in such turbulence is a goal for a future. 
In terms of star formation, we do not believe that turbulence 
is strongly imbalanced. In fact, in compressible media the 
imbalance decreases due to reflecting of waves from pre- 
existing density fluctuations and due to the development of 
parametric instabililites (see Del Zanna, Velli & Londrillo 
2001). 

12 Conclusions 

In this paper we consider the process of reconnection diffu- 
sion which arises in magnetized turbulence due to fast mag- 
netic reconnection. We consider GS95 model of strong tur- 
bulence and LV99 model of fast magnetic reconnection. We 
claim that in the presence of strong magnetic turbulence the 
notion of plasma frozen onto a magnetic flux tube becomes 
meaningless. Due to reconnection, magnetic field lines con- 
stantly change their identity inducing intensive mixing of 
plasmas. The intensity of this process is determined by the 
intensity of turbulence and does not depend on microphysics 
of reconnection, e.g. it does not depend on the collisional or 
collisionless nature of the small scale reconnection events. 

While reconnection diffusion is important for many key 
astrophysical processes, e.g. heat transport in plasmas, gen- 
eration of magnetic field etc., this paper is focused on star 
formation, where the flux freezing idea is at the core of the 
existing paradigm. Therefore the concept of reconnection 
diffusion alters the paradigm. The new concept explains the 
existing observational data that contradicts the theory based 
on ambipolar diffusion. In particular it can explain data on 
the decorrelation of density and magnetic field in diffuse in- 
terstellar media, fast removal of magnetic field in molecu- 
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lar clouds, higher magnetization of cloud core compared to 
envelope, properties of circumstellar rotation disks, indepen- 
dence of star formation rate on metallicity, possibility of star 
formation in ULIRGs, selected empirical relations etc. It 
also provide predictions for future observations and helps to 
bridge the gap between numerical simulations and actual star 
formation in galaxies and early universe. More studies of the 
new approach to magnetic field diffusion in star formation is 
necessary. 
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